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I. INTRODUCTION 

The theory of many-fermion systems is of obvious 
importance because of its ubiquitous applications in 
condensed-matter physics. It is also a very hard prob- 
lem, due to the difficulty of dealing with a macroscopic 
number of electrons which, in many materials, interact 
strongly enough to preclude treating the interactions as 
weak perturbations of noninteracting electrons. Histori- 
cally, there have been two main approaches to the prob- 
lem. Landau's Fermi-liquid theoryjEEl although not per- 
turbative, is phenomenological and self-contained in na- 
ture, and often a more microscopic approach is desirable, 
especially for investigations of the limitation of Fermi- 
liquid theory. Many-body diagrammatic techniques)^ 
on the other hand, are limited by their inherently per- 
turbative approach. Furthermore, it is neither feasible 
nor desirable to study all aspects of a microscopic model, 
except in some special cases. For many important appli- 
cations it suffices to keep only the low-lying excitations 
or soft modes of the system, which govern the behav- 
ior at long times and large distances. The most obvi- 
ous examples of such applications are phase transitions, 
but there also are many qualitative properties of entire 
phases that depend on the soft modes only. These observ- 
ables also have the attractive property of being universal 
in the sense that they depend only on basic symmetry 
properties of the system, and not on any microscopic 
details. Consider, for instance, a Heisenberg ferromag- 
net. The critical behavior at the Curie point, i.e., the 
phase transition from the paramagnetic to the ferromag- 
netic phase, is entirely governed by soft modes,!^ but so 
is the divergence of the longitudinal susceptibility every- 



where in the ordered phase.'^To describe and understand 
either phenomenon it therefore suffices to consider an ef- 
fective theory that takes into account only the relevant 
soft modes. In the case of the Heisenberg ferromagnet at 
temperatures T > 0, either a classical 0(3) (^'^-theory or 
an 0(3) nonlinear sigma model provides a suitable effec- 
tive theory.'"^ Such a theory can sometimes be constructed 
based on symmetry considerations only, but it often is de- 
sirable, and leads to greater physical insight, to start with 
a microscopic description and derive the effective theory 
by integrating out all non-essential degrees of freedom 
in some simple approximation that respects the crucial 
symmetries. The resulting "hydrodynamic" , i.e., long- 
wavclcngth and low-frequency, description of any system 
of interest is complementary to, and serves a very differ- 
ent purpose than, for instance, a first-principles approach 
that aims at a description at atomic length and energy 
scales. In order to obtain a well-behaved effective theory, 
care must be taken to keep all soft modes of interest and 
all other soft modes that couple to them. Integrating 
out soft modes leads to long-ranged effective interactions 
between the remaining degrees of freedom, and hence to 
non-local effective field theories that can be hard to han- 
dle technically. 

For noninteracting disordered electrons such an 
effective-theory approach was pioneered by Wegner.^^ 
Guided by an analogy with classical ferromagnets he 
constructed a matrix nonlinear sigma model capable of 
describing the instability of a (disordered) Fermi-liquid 
phase against the formation of an Anderson insulator, 
and allowed for the determination of the critical be- 
havior at the Anderson transition in terms of an e- 
expansion about the lower critical dimension for this 
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problem, d~ = 2. Wegner's original formulation was ex- 
panded upon by various authorsf^^ and eventually his 
sigma model was generalized to interacting disordered 
electrons fl2l This theory and its generalizations have been 
extensively used to describe the Anderson-Mott metal- 
insulator transitionli^^ These theories are fundamen- 
tally based on diffusive soft modes, and do not allow for 
the clean limit to be taken. This led to the strange sit- 
uation where a successful effective soft-mode theory was 
available for disordered electrons, but not for clean ones, 
although one might think that the latter poses a sim- 
pler problem.i^ This last conclusion is fallacious, how- 
ever, mostly because the number of soft modes is much 
larger in clean electron systems than in disordered ones. 
This was noted in Ref. '16' whose main objective was to 
provide a thorough derivation of Finkelstein's generalized 
nonlinear sigma model for disordered systems. While this 
theory did allow for the clean limit to be taken, the re- 
sulting field theory had non-local vertices and did not 
offer any obvious way to obtain a description as useful 
as the nonlinear sigma model for the disordered case. 
Another effective field theory for clean electrons was de- 
veloped in Ref. |T71 This theory is technically based on 
a supersymmetric matrix formulation and aims to in- 
corporate and unify various other approaches that dis- 
cuss the validity and limitations of Landau Fermi-liquid 
theory.liSHllI Jt introduces bosonic degrees of freedom and 
integrates out the fermions, as did Ref. (THl in contrast to 
Shankar's derivation of Fermi-liquid theory that applied 
renormalization-group (RG) techniques to a fermionic 
field theory.Eil The theory of Ref. [T7] has recently been 
refinecP^ and formulated both in a version that is for- 
mally exact and may be better suitable for numerical 
work than standard Monte-Carlo techniques and as 
an effective field theory for low-lying excitations .l^^ The 
latter was used for an analysis of the low-temperature 
specific heat in two-dimensional systems. 

In the present paper we take a different approach to 
this problem. It is similar in spirit to Refs. 1171271 but 
rather different technically. Our physical motivation is 
the desirability of a broadly useful effective theory that 
is capable of describing widely different phenomena such 
as ferromagne tic or dering in non-s-wave angular momen- 
tum channels ,'2Sl22] quantum ferrimagnetic order, which 
has received little attention so far, a possible breakdown 
of Fermi-liquid theory due to a vanishing quasi-particle 
weight,'^ and, in addition, novel universal aspects of 
the Fermi-liquid phase.'sSl We are interested in a local 
field theory that is not perturbative with respect to the 
electron-electron interaction, but rather can be controlled 
by means of a systematic loop expansion. This will pro- 
vide a basis for an application of RG techniques that 
allow for a resummation of the expansion in powers of 
the loop parameter, and thus allow to go beyond per- 
turbation theory in a controlled fashion. This not pos- 
sible with in the framework of traditional perturbation 
theory,!^ which provides no information about the struc- 
ture of the renormalized theory. Technically, the theory 



we will develop is based on a generalization of Ref. [121 
The reason for the non-locality of the latter was its re- 
striction to density or s-wave {£ — 0) modes, with soft 
modes in higher angular-momentum channels being effec- 
tively integrated out. Formulating the theory in terms of 
phase-space variables allows one to keep the soft modes 
in all angular-momentum channels. This leads both to a 
local theory and allows for natural description of order 
in ^ 7^ channels. Integrating out the massive modes 
in a tree approximation leads to an effective action that 
explicitly keeps all of the soft modes. Its structure is dif- 
ferent from that of the generalized nonlinear sigma model 
for the disordered case. It can be constructed to any de- 
sired order in the soft degrees of freedom, which in turn 
allows to perform a loop expansion to any desired order. 

The organization of this paper is as follows. In Sec. [ll] 
we formulate our model in terms of fermionic fields. We 
then bosonize the theory by constraining bilinear prod- 
ucts of ferniion fields to bosonic matrix fields Q and in- 
tegrating out the fermions. We then discuss symmetry 
properties of the Q-fields and representations of observ- 
ables in terms of Q-correlation functions. In Sec. |III| we 
identify the soft modes of the action by means of a Ward 



identity. In Sec. |IV|we expand the Q-field theory about 
a saddle-point solution that describes a Fermi liquid, and 
separate the soft fluctuations from the massive ones by 
means of the Ward identity. The massive degrees of free- 
dom are then integrated out in a tree approximation to 
derive an effective theory entirely in terms of the soft 
modes. We also make contact with the density formula- 
tion of Ref. UM In Sec. El we discuss universal features of 
the energy-dependent density of states (DOS) in a Fermi- 
liquid as a simple application. In Sec. |VI| we summarize 
and discuss the effective theory. We finally discuss var- 
ious future applications of the theory. A discussion of 
the effects of quenched disorder, and a pedagogical dis- 
cussion of the 0{2) nonlinear sigma model that stresses 
analogies with the current theory, are relegated to two 
appendixes. 



II. FIELD-THEORETIC FORMULATION OF 
THE INTERACTING FERMION PROBLEM 

A. Fermionic formulation 

Our starting point is a description of interacting 
fermions in terms of Grassmann fields. The partition 
function can be writterPl' 



(2.1) 



Here 'ip and ip are Grassmann- valued fields , and 

is the Grassmannian integration measure. The action 5' 

has the form 

•' a 

(2.2a) 
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The first term describes noninteracting electrons with 
chemical potential /x. For simplicity and definiteness we 
will assume a single parabolic band, i.e., 



e(V) = -VV2too 



(2.2b) 



with nic the effective electron mass. If desired, the model 
can easily be generalized to include a nontrivial band 
structure instead. We use a (c?+ l)-vector space-time no- 
tation with X — {x, t), J dx — Jy dx j^^"^ dr. x denotes 
the position, r is the imaginary-time variable, V and T 
are the system volume and the temperature, respectively, 
and a =t,i= +, — is the spin label. We use units such 
that h = fee — 1- '5'int describes an electron-electron 
interaction via a two-body potential v{x), 



S\nt — 



y<-lpat{xi)i^a^{x2)i^a^{x2)i^ax{xi) . (2.2c) 



In a microscopic theory the potential v{x) would be 
the Coulomb interaction. Here we assume for simplic- 
ity that the theory has already been renormalized to 
take into account screening, so v(x) represents a stat- 
ically screened Coulomb interaction, or a some simi- 
lar short-ranged model interaction. Equation (2.2c) de- 



scribes the interactions of number-density fluctuations at 
all wavelengths, and in a microscopic description this 
is the only interaction term there is.'^^ Integrating out 
fluctuations at short wavelengths in order to generate 
an effective long-wavelength theory generates interaction 
amplitudes in the spin-triplet channel in addition to the 
spin-singict one, and in the particle-particle channel in 
addition to the particle-hole hole one, and all of these ef- 
fective interaction amplit udes appear in all angular mo- 
mentum channels .'^^'^IIM! xhe theory we will develop is 
general, and any desired interaction amplitudes can be 
kept. However, for the sake of transparency and simplic- 
ity of the formalism we will restrict ourselves to a finite 
number. Some of the most interesting applications of the 
theory are related to instabilities of the Fermi-liquid state 
due to a density-density interaction, and to magnetic or- 
der of the ferromagnetic or magnetic nematic variety. 
Furthermore, the soft modes in the particle-particle chan- 
nel become massive in the presence of a magnetic field, 
as they do in disordered systems,'^! and hence can be 
suppressed experimentally. Accordingly, in the remain- 
der of this paper we keep only interaction amplitudes in 
the particle-hole channel, and further restrict ourselves 
to the spin-singlet s-wave {(. = 0) and the spin-triplet 
s-wave {£ = 0) and p-wave (£ = 1) channels. Gener- 
alizations to other interaction channels, for instance, the 
spin-triplet d-wave channel that has been discussed in the 
literature,12s3 are straightforward. It needs to be stressed, 
however, that a restriction to low-angular-momentum 
channels represents an approximation over and above the 
low-energy restriction that is central to our approach, 
and that a complete low-energy effective theory needs to 
keep all angular momentum channels. 



We define Fourier transforms of the fermionic field 
EE 7/;„,(fe) - I dx e-'^'^^^ix), 



dx 



-Mx), (2.3) 



where k = {k,ujn) is a four- vector that comprises a 
wave vector k and a fermionic Matsubara frequency 
u!n = 27rT(n -I- 1/2), and kx = k ■ x — a;„r. With the 
restrictions explained above the action then reads 

S = '^^l:„{k) [iujn-k^/2m^ + lJi\ ipaik) + Sint, (2.4a) 



with 



if) T 
2 V 



.(0) 



T 



2 V 

2 V 



q i—1 

q 



Here q — (q, r2„) comprises a wave vector q and a bosonic 
Matsubara frequency ri„ = 2T:Tn. ri°\ Tf\ and F^^^ 
are the interaction amplitudes in the particle-hole spin- 
singlet s-wave, spin-triplet s-wave, and spin-triplet p- 
wave channels, respectively. They are related to the pa- 
rameters Fq, Fq, and Ff , respectively, in Landau Fermi- 
liquid theory.l^ denotes a sum over wave vectors that 
is restricted to |<7| < A with a cutoff wave number A. 
The long-wavelength properties we are interested in do 
not depend on A. The electron number density n, spin 
density rig, and spin current density jg" are given in Eqs. 
(2.5) below, and for completeness we also list the number 



current density j"^ 



n{q) 



EE^^-(^ + '?/2)^-(^-9/2), (2.5a) 

k a 

Y.^'-Y. ^-^^ + Mk - g/2), (2.5b) 

k <y 

Y Y '^'^i + (c^j)'Ti-T2 V'rTs {k - ql'2.), 



k 0'1,0'2 



(2.5c) 

E '/'<Tl(fc + 9/2)(CTi)<Ti<T2V'<T2(fc-'7/2). 

k cri;(72 

(2.5d) 



The <7i {i = 1,2,3) are the Pauli matrices, and the ka 
[a = 1,2,3) are the components of the wave vector k. 

T^s^ > describes a repulsive density-density interaction, 
and Fj'^ > describe ferromagnetic interactions in the 
s- and p-wave channels, respectively. 
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B. Bosonic formulation 



with an adjoin' 



1351 



1. Phase-space degrees of freedom 

Our goal is to rewrite the action in terms of bosonic 
matrix fields. For this purpose it is convenient to define 
a bispinoJli^ 



1 
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-V^nt(x) 



(2.6b) 



V2 



-Ipntix), -Ipniix), Ipniix), -i^ntix)) ■ 



Here C„ 



i{Ti (g) S2) ^nm is the charge conjugation 



matrix in the spin-quaternion space spanned by (g) sj 



= 0, 1, 2, 3) with Ti 



-iaj with ctq the 2x2 



(2.6a) unit matrix and cri^2,3 the Pauli matrices. Explicitly, 



To 

So 



1 

1 

1 
1 



Tl 



Si 



We now define a bilinear tensor product 



Bnmix.y) 



;{x) 



-i 
-i 

i 

1 



® rimiy) 
--0„f (a;)'0„^(y) 
-'0n;(a;)V?mt(y) 

'0n;(a?)V^mt(y) 

■■!/'nt(a;)V'mt(y) 



T2 



S2 



-1 

1 

1 

-1 



-i^n^{x)ljj_m^{y) 
-tni{x)tmiiy) 
ll)_ni{x)lpjni{y) 
-1pnt{x)lpml{y) 



T3 
S3 



-'4j„f{x)lpmliy) 
-tnlix)'ipmliy) 

ip_niix)ipmi(y) 

-1pnt(x)lpml{y) 



-i 

i 

1 

-i 



i^nt{x)-lpmt{y) \ 
'tnlix)'ipmt(y) 
-ltnl{x)lpmt{y) 
'0nt(a;)V'mt(y) / 



(2.7a) 
(2.7b) 



(2.8) 



and its Fourier transform 



Bn^{k,p)^^ I dxdy e-'^-^+'P y Bnra{x,y)- (2.9a) 

The 4x4 matrix _B„m(fe,p) can be expanded in the spin- 
quaternion basis defined above, 

3 

B^^(k,p) = lB,^^{k,p) (r, ® s,) . (2.9b) 



i,r—0 

It is further useful to define 



B„ 



(fc;g) 



B„^{k + q/2,k-q/2) 



(2.9c) 



with analogous definitions for other objects that depend 
on two wavevectors. All bilinear products of the fermion 
fields Tp and ip can be written in terms of B, and in par- 
ticular all terms in the interacting part of the action, Eq. 
(2.4b), can be written in terms of products of the B. An 



i nsp ection shows that in the spin-quaternion basis, Eq. 



(2.7), the matrix elements ^B and ^=1.2-^ describe the 



particle-hole and particle-particle channels, respectively, 
and the '^^i? describe the spin-singlet channel while the 
describe the spin-triplet channel. Taking mo- 
ments of Bnm{k; q) with respect to k generates variables 
in different angular momentum channels. For instance, 



k 



ik;q) 



(2.10a) 



defines s-wave or density degrees of freedom, 

B^^l{q)^Y.^B^^{k-q) (2.10b) 



defines current degrees of freedom, etc. 

We note that the introduction of bispinors, and of ma- 
trices with spin-quaternion valued matrix elements, is 
necessary in order to handle both particle-particle and 
particle-hole degrees of freedom within the same frame- 
work, as it allows to form bilinear products ipip and V'V' in 
addition to i/j'f/', see Eq. ( |2.8| ). For purposes that involve 
the particle-hole channel only one can restrict oneself to 
ordinary spinors, in which case the matrix B^ Eq. (2.8 1, 
becomes a 2 x 2 matrix. Another observation is that 



B constitutes an overcomplete representation of bilinear 
fermion degrees of freedom since all of the matrix ele- 
ments of B are not independent. We will come back to 



this point in Sec. IIB 3 below. 



2. Q-matrix field theory 

Our next step is to constrain the matrices B in the 
interaction terms to a classical matrix field Q by means 
of a Lagrange multiplier field A. The fermion fields then 
enter the action only bilinearly and can be integrated out 
excatly. This way we obtain an effective action A that 
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depends on Q and A according to 

Z ^ j D[4>,i;] e^['^''^I ^ j D[n\ e^l"! 

D[Q,A] e^l'^.A] . (2.11) 



to the phase space variables, Eqs. (2.8 2.91 yields the 



Here and it what follows Tr denotes a trace over all de- 
grees of freedom, including the continuous position in real 
space, while by tr we will denote a trace over all discrete 
degrees of freedom that are not explicitly shown. 

In order to describe phenomena in angular-momentum 
channels higher than the s-wave one, or even in order 
to keep non-s-wave interaction constants, one obviously 
needs to apply this pr ocedure to the phase-space vari- 
ables defined in Eqs. (2.8 2.9). As we will see later. 



this leads to a local field theory that allows for a system- 
atic loop expansion about a saddle-point solution that 
describes a Fermi liquid. For certain purposes that re- 
quire a density-channel interaction only it is technically 
advantageous to formulate the theory in terms of density 
variables, Eq. (2. 10a I, even though this leads to a non- 



local theory, i.e., one where the vertices diverge in the 
limit of small wave numbers. The starting point for such 
a formulation was given in Ref. 1161 and we will revisit 
and further develop this method in Sec. |IV E| 



Application of the procedure sketched in Eq. (2.11) where 



following formal expression for the effective action: 



A[QM = ^0 + Tr A^Q + A„t[g] 



Here 



where 



A= ^TrlnG- 



G-^ = Gn ' - iA 



(2.12a) 



(2.12b) 



(2.12c) 



is the inverse Green operator, with 



X [iujn + V^/2mc + ^i] (2.12d) 

the bare Green operator. If a nontrivial band structure 
is desired, — V^/2TOe should be replaced by an appro- 
priate energy function e(V), see Eqs. (2.2). Rewriting 
the interaction part of the action, Eq. (2.4b), in terms of 



the and constraining the latter to Q by means of the 
functional delta-constraint we find 



(0,s) 



A- t — 



- A^ 



(0,t) 



-A' 



(i,t) 



(2.13a) 



A- 



(0,s) 



A^ 



A^ 



TF 



(0) 



2V 



r,i5=0,3 "l."2 

"3, "4 



k,p q 



(0,t) _ 



TV 



(0) 



t_ 

2V 



X! X! '^^Vi-n2.«4-n3X!X! *''('^'^'-®*^)'5«in2(fc;q))(Tr«'Si)(9n3„4(p;-q) 

fc,p q 



r,s=0.3 1=1 "i."2 

"3,>M 



(2.13b) 



(2.13c) 



(i,t) _ 



TF 



(1) 



^ 

2V 



^ ^ ( ) ^ ^ ^ ^ '^ni-n2,ri4-ri3 
r,s=0.3 1=1 "i."2 

"3,>M 



^ fc • p ^ tr (^(r^ (g) Sj) Qmn2 [k; q)^ (r^ ® Si) Q^^m (p; ~q) ^ 
q 

(2.13d) 

I 



3. Symmetry properties, and representation of observables 



We now derive some useful symmetry properties of the 



Q-matrices. B as defined in Eq. (2.81 is self-adjoint un- 



der the operation defined in Eq. (2.6b I. Q inherits this 
property, so we have 



Q+ = C' 



c 



(2.14a) 



In the spin-quaternion basis defined in Eq. (2.7) this im- 
plies 



+ 
+ 
+ 



+ 
+ 
+ 



+ 
+ 
+ 



(2.14b) 



(2.14c) 



(2.14d) 
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Here the symbols 



etc. denote a factor of +1 



for r = and a factor of —1 for r — 1,2,3, and anal- 
ogously for i. We have made use of these relations in 
order to write the interacting part of the action in the 
form of Eqs. (2.13). They imply that all of the Q- matrix 



elements are not independent. In a model with N Mat- 
subara frequencies, only iV(iV-|-l)/2 matrix elements are 
independent. We will later choose these to be the ones 
with n > m. 

In terms of the Q-matrices the observables listed in 
Eqs. (2.5) are given by 



m r=0,3 k 

(2.15a) 

j„(q,zao = E E E^"°Q™.„+Jfe;g) , 



7n r— 0,3 



(2.15b) 



EEH^^^'E°^?"..™+«(fe;'?)' 

k 

(2.15c) 



m r— 0,3 



k 

(2.15d) 



m r— 0,3 



These are just examples of spin-singlet and spin-triplet 
observables in the £ = and 1 = 1 channels, respec- 
tively. Clearly, and desired observables can be expressed 
in terms of the Q. 



4- Correlation functions 

Physical correlation functions can be written in terms 
of Q-correlation functions by keeping appropriate source 
terms in the action while performing the transformation 
to the Q-matrix variables. For instance, for the Green 
function Gnaix - y) = {^pnaix) tpnaiv)) we obtain 



Gna{x-y) =2*^ {{To+iT3)(E){soTcris3)) {Qnn{x,y)) 

(2.16) 

with a = "1,1 = +1,-1. This can also be read off of 
Eq. (2.8) directly by keeping in mind the isomorphism 



between B and Q. For the density of states as a func- 
tion of the distance lj from the Fermi surface this implies 



7V(^) = -Re(gg„„(a;,a^)) 



(2.17) 



where we have used the symmetry properties expressed in 
Eqs. (2.14 1. Similarly, the number density susceptibility 



X, and the spin density susceptibility tensor xl'' , can be 
written 

X{q,^n.n) = ler^ ^ ^(^(AQ)™,_„,™,(fc;q) 

mi, 7712 r— 0,3 fc,p 

x;(AQ)„,,„3+„(p;-q)) , (2.18a) 
Xng^i^n) = 16T^ ^ ^(;(Ag)„,,_„,„,,(fc;q) 

mi ,m2 r=0,3 k,p 

xJ,(AQ)™,,„,+„(p;-q.)) , (2.18b) 



Here AQ = Q - (Q), and, in Eq. ( pT8b| , ij = 1,2,3. 
Other correlation functions can be expressed analogously. 



III. IDENTIFICATION OF SOFT MODES: 
BROKEN SYMMETRY, AND A WARD 
IDENTITY 

We are interested in separating the degrees of free- 
dom represented by the Q-matrices into soft and massive 
modes. To identify the soft modes we derive a Ward 
identity that relates two-point Q-correlation functions to 
other quantities. We first do so for noninteracting elec- 
trons and then discuss the effects of interactions. 



A. A Ward identity for noninteracting electrons 

Let us consider transformations of the bispinors de- 
fined in Eqs. (2.6), 



Vn{x) - 



dy {x,y)iim{y) , 



(3.1a) 



where the operator T'^^ defines non-local infinitesimal 
rotations in frequency space, 

TitHx,y) = iTo®So)tt2ix,y) , (3.1b) 
^V±{x,y)+0{ip') . (3.1c) 



Here the 



V±{x,y) 



{x,y) ± (l){y,x)] 



±^±{y,x) 



(3.1d) 



are even and odd combinations, respectively, of a non- 
local rotation angle <j). The matrices obey (T^^')-^ = 
(f(±))-i, and (r(±))'^Cr(±) = C wi th C the charge 
conjugation matrix defined below Eq. (2.6b I. For fixed 



TLi and n2 they form an SO (2) subgroup of a much larger 
symplectic group (Sp(8iV, C) for a system with 2N fre- 
quency labels; N positive ones, including zero, and N 
negative ones) , see Ref. TSl. The Q-matrices transform as 



Q(a:,y)^(T(±)g(r(±))^ 



(3.2) 
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Explicitly we find with 

Qnm{x,y) Qnrn{x,y) + 5Qnrn{x,y) (3.3a) 



5Qnm{x,y) j dz |(/7±(y,2;) [5,nm Qnn^ix, z) ^ S^n^ Q nm{x , z)]+Lp ±{x , z) [S^m Qn.27n{z,y)TSnn2 Qnim{z,y)]^ , 

(3.3b) 

and in particular 



{x,z)(p±{z,y) 



(3.3c) 



The Lagrange multiplier field A transforms as Q does, 
on account of the bilinear coupling between the two. Of 
the three terms in the action, Eqs. (2.121, the second 



one is invariant under these transformations, but stnd 
Aint are not. Focusing on noninteracting electrons for 
the time being, we find Aq — > Aq + SAq with 

6Ao = ^ [iri„i_„2 + • q/'^c] tr G„2„i (fc; q) 

k,q 

xip±{k-q/2,k + q/2) . (3.4) 

Here G is the inverse of the Green operator defined in 
Eq. (|2.12c|). 



A Ward identity can now be derived by standard 
techniques.^ We introduce a matrix source field J for the 
Q, consider the generating functional 



Z[J] = / D[Q,A] 



„^o+Tr(A^Q)+Tr (JQ) 



(3.5) 



perform the infinitesimal rotation defined by ip±, differ- 
entiate with respect to J, and put J = 0. This way we 
obtain a Ward identity 



{SAo (9„i«2 (a:, y))Ao + {^Qnin2 {x, y))^^ = . (3.6) 



From Eqs. (3.4) and (3.6 1 we see that this relates cor- 



relation functions of the structure (trGQ) to (Q). The 
former can be rewritten in terms of {Q Q) by generaliz- 
ing the generating functional given in Eq. ( |3.5[ ). Since 
the Q are isomorphic to B, Eq. (2.8), we can write the 



source term JQ — xJQ -I- (1 — x)JB with an arbitrary 
real number x. The generating functional then becomes 



Z[J] = / D[Q,K] e5 



1 Tr ln[G"i+i(l-a:)J^]+a;Tr(jg)+Tr(A^Q) 



(3.7) 



Note that this is independent of x, and that by choosing x = 1 we recover Eq. (3.5). By differentiating with respect 

(3.8a) 



to J, choosing x = Q and x = \, respectively, and putting J = we obtain an identity 

{Gn2nAx2,Xi)) = -2i {Qn^n2{xi , X2)) ■ 

Differentiating twice with respect to J we find 

{Gn2ni{x2,Xi) QnanA^3, X4)) = -2i {Qnin2{xi, X2) QnsnA^d, X^)) 



(3.8b) 



Adding them yields the Ward identity in its final form 



We now differentiate Eq. (3.6) with respect to (f> in Eq. (3. Id). This yields two identities, one for ip^ and one for ip_ 

(S]Q„,„,(fc-g/2))-(SQ„,„,(fc + 9/2)> 



(3.9) 
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B. Soft modes in noninteracting electron systems 



Let us discuss the Ward identity for noninteracting electrons, Eq. (3.9). By using Eq. (3.8a I we see that the 
right-hand side of Eq. (3.9 1 is determined by 



2i {Q^^ik ± q/2)) = (G„„(fc + q/2)) 



1 



— £,k±q/2 



q/2) 



l/^fc±<j/2 



(3.10) 



Here — /2mp — /i, and the last term on indicates 
the limiting value of the Green function as aj„ approaches 
zero from above or below, respectively. We see that if ni 
and 71,2 have opposite signs, then the numerator on the 
right-hand side of Eq. (3.9) goes to zero as fli-2 — 5> and 
q — >■ 0, as docs the denominator, whereas the numerator 
remains finite if ni and 712 have opposite signs. In the 
latter case, the leading term for small and small q 

is 



(oQ„,„.(fc;g)oQ„,„.(p;-'z)) 

J_ 27ri^(gfc)sgn r^ni-«2 
16 ifl„^-n2 + k ■ q/rric 



+ 0(1) . (3.11) 



We see that there is an infinite number of soft modes 
that can be obtained by taking all possible moments of 
Eq. (3.11 ) with respect to the center-of-mass wave vector 
k. This identifies the matrix elements Qnm with nm < 
as soft degrees of freedom, whereas the Qnm with nm > 
are massive. Physically, the soft modes are particle- 
hole excitations with a linear frequency-momentum re- 
lation. In the absence of interactions they are massless 
at nonzero temperature as well as at T = 0. 



Equation (3.11) is a generalization of the Ward iden- 



tity derived in Ref. 1161 which corresponded to the zeroth 
moment at q = 0, 



TriVp 



q=0 



(3.12) 



Here Q'-°^( g) — J2k Qi^', q) is defined in analogy to 

Eq. ( |2.10a[ ), and Np is the free-electro n den sity of states 

per spin at the Fermi level. Equation (3.12) is the clean 



analog of an identit y that has been discussed before for 
disordered electrons PISISlHalljn the presence of quenched 
disorder this is the only soft mode; all higher moments 
acquire a mass, see AppendixjXj This is the reason why it 
is much harder to construct a soft-mode theory for clean 
electrons than for disordered ones: In the former case 
there are many more soft modes. 

The Ward identity shows only that the matrix elements 
oQnm with nm < are soft. However, it is easy to see 
that the same statements holds for the IQ with arbi- 
trary values of i and r. The reason is that the (IQIQ) 
correlation functions are related to (qQ qQ) by means of 
symmetries that are not broken. This argument is iden- 
tical to the one give in Ref. 16 and we do not repeat it 
here. 



C. Interaction effects 

We now need to address the question of what hap- 
pens if the electron-electron interactions are taken into 
account. To this end we need to add the variation of the 
interacting part of the action, SA-mt, under the trans- 
formation given in Eqs. (3.1), to Eq. (3.6). 6Aint is 



quadratic in Q, and the Ward identity now relates the 
2-point function on the left-hand side of Eq. (3.9) to a 



1-point function and a 3-point function. The technical 
details are cumbersome and have been derived in Ref. 
[T6l but to sh ow t hat the particle-hole excitations iden- 
tified in Eq. (3.9) remain soft it suffices to discuss the 



structure of the generalized Ward identity, as was shown 
in Ref. 39 , In what follows we recapitulate and adapt the 



argument give there. Equation (3.9) gets generalized to 



-8i (iri„j_„2 +k ■ q/ml) Dni„2,„3ni{k,p;q) 

(3.13a) 

Here m* is the effective electron mass that is renormal- 
ized by the Fermi- liquid parameter , and 



Nn^n.^{k,q) = 7rZfc(5(^fc)sgnfJ„j_, 



(3.13b) 



with Zk the quasi-particle weight that rcfiects the change 
from free-electron Green functions to physical ones. W is 
the contribution from SA^t ■ General considerations show 
that there are two contributions to D that are character- 
ized by different frequency structures. 



D 



[k,p; q) 



^riins ^n^rii ^i'Jn2 ^iV'i ?) 



Dnin2,n3niik,p; q) , 

(3.14) 



and analogously for W. Here the superscripts (dc) 
and (c), refer to the disconnected and connected con- 
tributions, respectively, in a representation in terms of 
fermionic diagrams. These contributions must respect 
the Ward identity separately, and either one diverging 
in the limit of small frequencies and waves numbers will 
ensure the existence of soft modes. It therefore suffices 
to discuss D'^'^'^\ The Ward identity for the latter reads 

- 8i (in^.-n^ + k ■ q/ml) D^^Xi^^P'^ ^ 

4,piV„,„,(fc,q) - Wifl{k,p-q) . (3.15) 
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Now analytically continue to real frequencies G M, 
iflni-n2 — + iO, and consider the limit fl 0, q ^ 0. 
The only way for D^"^^) at t = to remain finite in 
this limit is for N and M^^"^^) at fl = q = to cancel. 
However this is in general not possible: A'^ approaches a 
nonzero limit for — > 0, g — > 0, and vanishing interac- 
tions, whereas W vanishes for noninteracting electrons. 
These two objects are thus different functions of the in- 
teraction, and their values at zero frequency and wave 
number can vanish at most for special values of the in- 
teraction. Notice that no such a cancellation is necessary 
to give a mass at nonzero temperature. For sim- 

plicity, consider the zeroth moment with respect to k and 
p only, and put q = 0. Then we have 



can be seen from Eq. ( 2.18a). Rather, the infinitely many 



SiflD^'^''^ =N- M^^'*") 



Now let 



£)(dc) 



(3.16) 



(3.17) 



with 1/t0 a phase relaxation rate. As long as fir^ 3> 1 
for T ^ and ^ 0, the right-hand side of Eq. ( [sle] ) 
can be expanded in powers of l/fir^, with the leading 
contribution being a constant. As long as l/r^ — for 
T — >■ such a regime always exists. 

We conclude that the soft modes identified in Sec. 
|III A[ namely, the Qnm with nm < 0, remain soft in the 
presence of interactions. The Ward identity shows this 
for the matrix elements [jQ in the spin-quaternion basis, 
but additional symmetries of the action ensure that it 
is true for all of the \.Q^^ This is consistent with both 
perturbation theory and with Fermi-liquid theory, which 
ensures that the number and nature of the soft modes 
cannot change as long as the system remains in a Fermi- 
liquid phase. The symmetry is broken, and the Gold- 
stone modes exist, as long as there is a non-vanishing 
quasi-particle weight Zp.. The only way the soft-mode 
spectrum can change is for the quasi-particle weight to 
vanish, which is to say that the system enters a non- 
Fermi-liquid state. We note that these considerations 
do not preclude using the soft-mode theory to study a 
breakdown of the Fermi liquid, which is signalized by a 
vanishing quasi-particle weight. They just indicate that a 
non-Fernii-liquid phase, where the symmetry is restored, 
has different properties and requires a different effective 
theory than the Fermi-liquid phase, but the theory we 
are about to construct will be valid everywhere in the 
latter, up to any possible transition where the quasi- 
particle weight vanishes. This is important for certain 
applications of the theory.'^Sl 

We finally note that while the right-hand side of Eq. 
(3.11) resembles the structure of the Lindhard function. 



the soft modes are not related to particle-number con- 
servation, and are not restricted to the density channel. 
Indeed, the density susceptibility cannot even be con- 
structed from the left-hand side, as doing so would re- 
quire three independent frequencies, rather than two, as 



soft modes are the Goldstone modes of a spontaneously 
broken continuous symmetry, namely, the rotations in 
frequency space expressed by Eqs. (3.1). 



IV. EFFECTIVE SOFT-MODE THEORY 



A. Saddle-point solutions 



We now return to the action, Eqs. (2.12 2.13), and 



consider saddle-point solutions. There are many saddle- 
point solutions that have different symmetry properties 
and describe different physical states. Here we give just 
two examples. The first one is a saddle point that de- 
scribes a Fermi liquid, and the second one reproduces the 
Stoner theory of ferromagnetism. Saddle points that de- 
scribe superconducting states, magnetic states with non- 
s-wave order parameters, or, if a nontrivial band struc- 
ture is taken into account, antiferromagnetic states, can 
be constructed analogously. 



1. Fermi-liquid saddle point 

The saddle-point equations are obtained by minimizing 
the action with respect to A and Q. They read 



~^7m{x.y) 



6Qnni{x,y) 



A„t[Q] ■(4.1b) 



To find a saddle-point solution that describes a Fermi 
liquid we make an ansatz 

Qnmix^y) = {tq<» So)SnmQn{x - y) , 

A^^^ix, y) = (ro ® sq) (5„„ A„(a; - y) . (4.2) 



Performing the derivative in Eq. (4.1a) and going into 



Fourier space we find Qn and A„ as the solution of the 
equations 



Qn{k) 

A„(fe) 



1 



(4.3a) 



Qmip) 



(4.3b) 



We see that —2iQn{k) equals the saddle-point Green 
fimction, in agreement with Eq. (3.8a), with A the self 



energy. The latter represents the spin-singlet interaction 
ri°^ in Hartree-Fock approximation. The factor e*'^'"° 



in Eq. (4.3b) is the usual convergence factor that re- 



solves the ambiguity inherent in the equal-time Green 
function.^ This Hartree-Fock saddle-point solution is a 
generalization of the saddle-point solution considered in 
Ref. Uni to a non-local Green function. 
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2. Stoner saddle point 



Gaussian fluctuations 



To illustrate the existence of other saddle points, and 
the general versatility of the theory, we now consider a 
saddle-point ansatz 



Qwmi^^y) = ^12 {to® Sn)Gn{x- y) 
+ ® S3.) Fn{x - y) 



(4.4a) 



-(to so)«S„(a; - y) 
+{T3(g)S3)iA„{x-y) .(4.4b) 



The saddle-point equations (|4.1|) now yield 
G„(fc) 



Sn(fc) 
A„(fc) 



4 

A = 



(4.5a) 

[g+{k)-g-{k)] , (4.5b) 

4zr(")i5:r^e^-°G„.(p), 

p m 

(4.5c) 
'"""i^™(p) . 

(4.5d) 



-Air. 



(0) 



V 



Here 



1 



(4.6) 



li^n - Cfc ± A - S 

are Green functions whose self energy contributions de- 
scribe the particle-hole spin-triplet and spin-singlet in- 
teractions in Hartree-Fock approximation. This saddle- 
point solution is the clean limit of the on considered in 
Ref. |3ni generalized to the case of non-local Green func- 
tions. It is equivalent to Stoner theory, which can be seen 
as follows. Gn{k) and Fn)k) obey the equations 

(iWn - a - S) g„(fc) + A j:„(fc) - 1 , (4.7a) 
(«w„-a-S) J-„(fe)+Ag„(fc) - . (4.7b) 

If we absorb the constant spin-singlet self energy E into 
a redefinition of the chemical potential we obtain from 



Eqs. (4.7) the equation of state in a form that is familiar 
from Stoner theory,Sll 

1 



1 = -2r 



(0) 



(4.8) 



which allows for a nonzero magnetization (which is pro- 
portional to A) provided iVpFj^^ > 1. 



B. Expansion about the saddle point 

In the remainder of this paper we consider fluctua- 
tions about the Fermi-liquid saddle point derived in Sec. 
|IV A 1[ an application of the theory to magnetically or- 
dered states will be discussed elsewhere.!^ 



We first consider Gaussian fluctuations about the 
Fermi-liquid saddle point. If we write Q = Q^^ + S O and 
A = A'^P -|- (SA, an expansion of the action, Eqs. (2.12 



2.13), yields 



A[Q, A] = y^^P + + , (4.9) 

where A^^ is the saddle-point action, A'^^^ denotes the 
Gaussian fluctuations, and A^ denotes the contributions 
that are of cubic or higher order in the fluctuations. For 
the Gaussian part we find 

^(2) ^ 1 Tr (Gsp 6~A Gsp 6 A) + Tr (SA^SQ) + Ant [SQ] . 

(4.10) 

Here 



(^sp)„,„ (fe,p) = (to (g) So) Snrn ^fc,p G„(fc) (4.11a) 



with 



GJk) 



1 



i^^n - - iA 



(4.11b) 



where A is the solution of Eqs. (4.3), is the Green func- 



tion in saddle-point approximation. With our local in- 
teraction amplitude A amounts only to a constant shift 
of the chemical potential, so G„(fc) can be taken to be 
the free-electron Green function. 

We can decouple 6A and 5Q at the Gaussian level by 
defining a field SA by 

;(M)i2 = ^ 012 ;(<5A)i2 + l{SQ)i2 ■ (4.12a) 

Here we have defined {SA)i2 = (5A)„j„2(pi,p2), and 
analogously for other fields, as well as 



;(<5Q)i2 -;(<5Q)i2 



(4.12b) 



and 



= Gi G2 = G„i(pi) G„2(p2) 

('$12 if (ni + l/2)(n2 + l/2)>0 
~ 1^12 if (ni-hl/2)(n2 + l/2) <0 . 

(4.12c) 

The frequency-restricted objects <I'i2 and 1^912 will be use- 
ful later. We also define 



W)i2 = {5Q)n,n2{-Pl,-P2) 



(4.12d) 



and analogously for SA. In terms of 5 A and SQ the Gaus- 
sian action then takes the form 



^2) = -4^ 




12 — ; (-50)12 



1 

1 

P12 



12 — ;(<5A)i2 



(4.13) 
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Note that the Gaussian SA propagator equals minus the 
6Q propagator^^' for noninteracting electrons, 



((5Ai2 5A34)_4(2) = -{SQl2SQ34).(2) 



(4.14) 



with A'-q'' the first two terms in Eq. (AAO) or ( |4.13 ). This 
will be important for the structure of the effective soft- 
mode theory and the loop expansion. 



2. The fluctuation action in terms of independent variables 



To go beyond Gaussian fl uctua tions we expand the 
Tr In term in the action, Eq. (2.12), in powers of SA and 



express the result in terms of 6 A and SQ by means of Eq. 
(|4T2a|. We find 



= - V — Tr(ap5A)'" = V AA^' 



m— 3 



4i 



- ^ Gr' tr [(Mi2 - <5Qi2)(M23 - <5g23)(<5A3i - SQ^^)] 



1,2,3 



-2 ^ Gr'G^'G^'Gjitr {{SA^^- bQ^^WA^z- SQ^-iWA^A- SQ^aWA^i- SQ^i)\ 



(4.15) 



1,2,3,4 



Notice that the vertices are given in terms of products of inverse Green functions, which are well behaved in the limit 
of small frequencies of wave numbers. This is thus a local field theory, in contrast to the formulation in terms of 
density modes in Ref. [1^1 

We now recall two features of the theory. First, all of the matrix elements of Q and A are not independent. We 
can choose the Qnm with n > m as the independent ones, and express those with n < m in terms of them by means 
of Eqs. (2.14). Second, we recall from Sec. Ill that the Qnm are massive degrees of freedom if n and m have the same 



sign, and soft ones if n and m have opposite signs. We express these facts in our notation by choosing as independent 
matrix elements 



(<5A)i2 = 



Pi2 if ni >n2 A (ni + l/2)(n2 + l/2) >0 

912 if rii > A 712 < , 

'Ai2 if ni > na A (ni + l/2)(n2 + 1/2) > 0) 

Ai2 if ni > A 712 < . 



(4.16a) 



(4.16b) 



In what follows we will absorb the frequency restrictions into the fields, so writing, e.g., qyi implies that ni > and 
7i2 < 0, etc. For the Gaussian action we obtain 



.(0) 



y- 



t 



^ n n J- T *P ^ p* «\ r *A *A 

ryi2 1-^12 ^ -"IS r-' 12 r-* 12 r^l2 r^l2 ^ 12 r-''-12 r^^l2 



i6Tr^ 
V 



r=0,3 1.2 

3,4 



1-2.3-4 



0„ Ot J Op Opt ,Q Opt I Op ot 
rHl2 r'dSA ' ^12 12 34 r'il2 34 12 ry34 



(4.17) 



where I12 = 1 — (5„^„,/2. Here we keep only the particle-hole spin-singlet interaction; the other interaction channels 
provide analogous contributions to the Gaussian action. Notice that the massless and massive degrees of freedom are 
decoupled in the noninteracting part of the action, but are coupled by the interaction via terms that are bilinear in q 
and P. 

The contributions to can also be expressed in terms of P, q, A, and A. Simple combinatorics yield for the cubic 
term 



AA^^^ = -4i Tr [G-^ {P - A) G-^ (P + P+ - A - A+) G-^ (P - A)+] 
+Tr [G-i (P + F+ - A - A+) G'^ {i G-^+ + i+ G'^ i)] 



(4.18) 
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or, schematically, 

and for the quartic one 

Ay^f'^) = -2Tr 



A^(3) ^ Q~3 [(p _ + (p _ A) 



or, schematically. 



AG-^{P-A)G-^{P- A) G-^ [P -A)G-^{P- A)+ 

+4 (P - A)+ (P - A)+ G-^ {P - A)+ G-^ {P - A) 
+3 G-^ {P - A) G-^ {P - A) (P - A)+ G^^ (P - A)+ 

+3 G-^ (P - A) G-^ (P - A)+ G-^ (P - A) G^^ (P - A)+ 

G"^(P + P+ - A - A+) G-^ (P + P+ - A - A+) G"^ G"^ g 

+G-1 (P + P+ - A - A+) G-^4G-^ (P + P+ - A - A+) G^^ 
-4Tr [G-i^G-i^G-i^G-i^+j , 



-8Tr 



A^(*) cx G~4 [(P - A)4 + (P - A) 



2 ^2 



G-1, 



(4.18') 



(4.19) 
(4.19') 



Here we have defined 



(4.20) 



and G^^P etc. denote matrix products with the ma- 
trix elements of G given by G12 = (512 Gi (tq (8) sq). 
The matrix elements of P are l.P„^j^^{pi,P2), see the 
left-hand side of Eq. (2.14c), and those of P+ are 
/ + \ / + \ 



^^n2ni(^P2,-Pi), see the right-hand 



side of Eq. (2.14c). The matrix elements of A, A+, q, 



q^, and A, A"*", respectively, are given by analogous ex- 
pressions. This implies for the expansion of the various 
matrices in the spin-quaternion basis 



r=0 1=0 
3 3 



Sj) , 

5 4) 



(4.21a) 
(4.21b) 



r=0 1=0 



where 



in the particle-hole channel, and 



r-^ 12 



r riin2 



(-P1 



-P2j 



(r = 0,3) (4.21c) 



(r = l,2) (4.21d) 



in the particle-particle channel, and are the her- 
mitian conjugates of and Si, respectively. Analogous 
relations hold for A, q, and A and their adjoints. In all 
cases the frequency ordering restrictions explained after 
Eqs. (|4.16|) apply. 



Expressions for the higher-order terms in IS.A in terms 
of P, g, A, and A can be constructed analogously. For 
instance, /S.A'^^'^ contains terms that have the schematic 
structures 



Ay^(5 



oc G-" [{P-Kf + {P-KYi' + {P-K)f] . (4.22) 



C. Integrating out the massive modes 



We now recall that P and A represent massive fluctu- 
ations, whereas q and A represent soft ones, and our goal 
is to construct an effective theory in terms of the latter. 
We cannot simply neglect P and A since they couple to 
q and A; rather, we need to integrate them out in an ap- 
proximation that respects the Ward identity discussed in 
Sec. |III| We let ourselves be guided in this process by an 
analysis of the 0(2) 0^-theory in Appendix [b| The proce- 
dure used below for our matrix theory is analogous to the 
one explained for the 0(2) model in the context of Eqs. 
(Bll ). To avoid unnecessary notational complexity, from 



here on we will restrict ourselves to a model with only 
one interaction amplitude, Fs*'^ = F, in the particle-hole 
spin-singlet channel. The spin-triplet (i = 1,2,3) and 
particle-particle (r = 1,2) degrees of freedom then just 
represent noninteracting electrons, and in the remainder 
of this paper we will neglect them. That is, we keep 
only matrix elements with r = 0, 3 and 1 = 0. It will be 
obvious how to generalize the theory to the presence of 
interactions in other channels.'^ 
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Screening of the interaction 



Equation (4.23a) thus simplifies to 



It is advantageous to first eliminate the coupling be- 
tween P and q at the Gaussian level. We see from Eq. 
(4.17) that this can be achieved by a simple shift of P. 



We define 



Up — " p 
12 ~ 12 



2TT 



J2 ^^r2!34^3-4,5-6°g56 • (4-23a) 



3,4,5,6 



Here 



M- 



12,34 



^13'^24<i'l2 

+<5l-2,3-4 



2(r/y)r$i2$34 



l-2TriE5,6'^l-2,5-6<>56 

(4.23b) 



is the inverse of 



Mi2,34 = '5l3'524'&i2^ 



2Tr 



Si- 



2,3-4 ■ 



(4.23c) 



Due to the frequency restrictions inherent in $12, see 
Eq. ( 4.12c| , the frequency sum in the denominator of 
Eq. (4.23b|) has no hydrodynamic content and can be re- 



placed by minus the static electron density susceptibility 
in Hartree-Fock approximation, Xst which in turn we can 
replace by Np: 



T 
V 



(^1-2,5-6 ^5 



Xst 



(4.23d) 



= 5'Fi2-2T7ci>i2:^^(5i-2,3-4j!g34 (4-24a) 

3,4 



Up 

12 



°(*7g) 



12 



with 



7 = r/(i-h2A^Fr) 



(4.24b) 



(4.24c) 



the statically screened interaction amplitude and opera- 
tors 

712,34 = (ro®so)2(T/y)7(5i_2,3-4 , (4.24d) 
^12,34 = (to «) So) (^13 (524 $12 • (4.24e) 

Analogously we define F and (p. These operators are 
all self-adjoint with respect to the matrix adjoint, i.e., 

In addition to decoupling P and q in the Gaussian 
action, this procedure introduces a new t erm q uadratic in 
q that combines with the q^ term in Eq. (4.17) to change 



the interaction amplitude F to 7. That is, the bilinear 
coupling between the massive modes and the soft ones 
leads to the screening of the interaction. The Gaussian 
action in terms of P, A, q, and A now reads 



5,6 



r=0,3 1.2 L V V^l 

3,4 

-8 E E^i2 [°Pi2 (sisS: 



<5l-2,3-4 ^, 

12 y 



^^^^ Oj A A ^ Ox Ox 

r934 ~ "I3O24 r^l2 r^. 

J fl2 



0\ 0\t 

34 



r=0,3 1,2 

3,4 



$12 



(^1-2, 3-4 



2ZI^ Opt 

V ) 



1 OA 0,1: 



34 - 013024 ^ ,.^M2 r^^Zi 
*12 



(4.25) 



The cubic and quartic parts, respectively, of the fluctua- 
tion action have the schematic form 



Ayt^^') oc G-^[[P - K + ^^qf 

+{P- K + ^-iq){q-\f] , (4.26a) 
A^(4) ^ [(P - A + $7(7)4 

+ {P- A + ^jqf{q~ Xf 

+{q - A)1 . (4.26b) 



More explicit expression are obtained by substituting 
Eqs. (|4.24|) into Eqs. (|4.18|) and (|4.19|. 



2. Integrating out the massive modes 

We now integrate out P and A in a saddle-point ap- 
proximation while keeping q and A fixed. This procedure 
is analogous to the one explained in Appendix [B] for de- 
riving the 0(2) nonlinear sigma model. In the current 
context it leads to an effective soft-mode theory that is 
related to, but has a different structure than, a sigma 
model. 

The saddle-point equations read 



= -16 /12 



l{^-'P)l,-'riml, +<5A^/<5°Pi 



12 ' 



= I6/12 "($"^A);2 -'5Ayl/50Pi2 



(4.27a) 
(4.27b) 



14 



Here we have used S/^A/SP =__5^A/5k, and VP is 
defined in analogy to 75 in Eq. (4.24b I. From Eqs. (4.27) 



we immediately obtain the useful identities 

(P-A)i2 = (<lfP)i2 

= (<l7A)i2 . 



(4.28) 



P is thus given in terms of A, which in turn can be ex- 
pressed in terms of q and A by solving Eq. (4.27b I. We 
will determine the action explici tly to (q'^), which re- 



quires A to 0(g ). From Eqs. (4.27b) and (4.26a) (or, 



more explicitly, from Eqs. (4.18) and ( 4.24[ ); note that 



only /S.A^^'^ contributes to A to this order) we find 



Al2 



-2i 



(($79) (<i>7<Z+))i2 

+ (($79+)G-i(<l79))i2 
+(($79) G-\^jq))i2 



+ 0{q'). 

(4.29) 



Terms of higher order in q can be obtained by means of 
an obvious iteration procedure. 

Fluctuations about this saddle point yield measure 
terms that reflect the Jacobian due to the change from in- 
tegrating over SQ and SA to obtain the partition function 
to integrating over q and A. For the 0(2) nonlinear sigma 
model this is demonstrated in Appendix [Bj While they 
are structurally important for preserving the symmetry 
of the action, these terms are not of practical importance 



for a loop expansion near the lower critical dimensionality 
d~ of the problem. The reason is that they lead to loop 
integrals that are less infrared divergent than those that 
result from the terms we have derived so far. In the 0{2) 
nonlinear sigma model they carry no gradients, see Eq. 
(B12), and hence are less infrared divergent near d^T = 2 



by two degrees than the terms contained in the nonlinear 
sigma model without the measure terms. The same state- 
ment is true for the generalized nonlinear sigma model 
that describes the Anderson-Mott metal-insulator transi- 
tion in the disordered interacting fermion problem, which 
also has rf- = For the current problem the lower 

critical dimension is d~ = 1, the dimension at and below 
which the Fermi liquid is not stable for any parameter 
values. This reduction of d~ compared to the disordered 
case is a consequence of the soft modes having a ballistic 
dispersion relation, as opposed to a diffusive one. For the 
same reason the measure terms will lead to loop integrals 
that, near d~ = 1, are less divergent by one power of mo- 
mentum than the ones we keep. For the sake of simplicity 
we will ignore these measure terms in what follows; it is 
clear how to derive them if that's desirable. 

D. Effective action 



Effective action to 0{q''') 



We are now in a position to write the effective action 
in terms of the soft modes q and A only. The relation 
between TP and 7A in Eq. (4.28) allows to write the 
Gaussian action, Eq. (14. 251), 



EE 

r-0,3 1,2 

3,4 



3"24 



"^12 



A 2r7\ t 

(51-2,3-4 I ,.^34 



1 



^13^24 r'^12 r'^34 

■^312 



"A 



2r7 

1^ 



^l-2.3-4°A 



34 



(4.30) 



This can be written more compactly 

_^(2) ^ -2Tr[q((^-i -7)g+] +2Tr (A(^-i A+) 

-2Tr (A7A+) . (4.31) 



Here A is given in terms of q and A by Eq. (4.29). Sim- 
ilarly, the other terms in the action, AA''"^\ AA'^ \ etc., 
can be expressed in terms of q and A by means of Eqs. 
( |4.24[ ), ( |4.28[ ), and ( |4.29[ ). The resulting action is the de- 
sired effective action in terms of the soft modes q and 
A. 

Note that this is a local field theory, i.e., all vertices are 
finite in the limit of small wave numbers and small fre- 
quencies. This can be seen most easily from Eqs. (4.26) 
and their generalizations: The terms with the highest 
power of Green functions in A^^"^ has the structure 
G-"$"(7g)" ~ G"(7g)". Since G ~ l/fi and 79 ~ Qq, 
this vertex is of 0(1), and all other vertices scale as a 



positive power of the frequency. 

One might consider analyzing this theory by means 
of a loop expansion in terms of both q and A. How- 
ever, this is not desirable. The reason is that the origin 
of A is a hard (i.e., delta- function) constraint that con- 
strains the q to soft bilinear fermion modes. Treating A 
perturbatively relaxes this hard constraint, which effec- 
tively introduces spurious soft modes. (We note that no 
such difficulty arises for A, which enforces a constraint 
for massive modes.) It therefore is desirable to treat A 
exactly. This can be achieved by eliminating A in fa- 
vor of diagram rules for an action that is formulated en- 
tirely in terms of q and ^. To see this let us consider the 
terms quadratic in q and A, i.e., the first two terms in Eq. 
(4.30) or (4.31). The Gaussian q and A-propagators are 
obtained by inverting the corresponding vertices. We see 
that the A-propagator is minus the noninteracting part 
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of the q-propagator, which is given by </3, and that q and 
A are not coupled in the Gaussian action. This implies 
that whenever ^ is contracted with ^ , the resulting prop- 
agator is effectively equal to the interacting part of the 
g-propagator, whereas contracting ^ with A results in the 
full g-propagator, as does contracting q with q^ A can 
thus formally be integrated out, despite the fact that it is 
a soft fluctuation, without causing undesirable features 



of the effective theory, such as nonlocality. The reason is 
that the only effect of A is to cancel well-defined contri- 
butions from other soft fluctuations. We thus arrive at 
the following effective soft-mode action: 

Aff = 4ff + A^fff^ + a4^3.) + 0(g5) . (4.32a) 
with a Gaussian part 



(2) 
eft 



~2Tr [q{r^-^)q+] 

^ XI $I°9i2 ('^13'524 - (51-2,3-4 °qli 



r=0,3 1.2 

3,4 



^12 



(4.32b) 



and non-Ga ussian parts Ay^[,^'' that each contain terms of 0(g"). The first two we obtain from Eqs. (4.18), (4.191, 
and (|4.24b|) as follows,E3 



A^, 



(3) 
cff 



-4iTr 
-4iTr 



G"i($7g) G"i ($7g + l>7g+)G"i ('^7 9+) 



(4.32c) 



A4^^)=-4Tr iG-' i+ iG-^+] 



-4i Tr 



G-i ($7A + <l7A+) G-^ a G-^ j 



+ G-1, 



-2 Tr (A7A+) 
-8Tr 



-4iTr 



-2Tr 



G-i($7(7 + $7(7+) G-i ($79 + $75+) G^^ (^G^^ ^+ + ^+ G'^ ^) 

+G"i ($7g + ^^q+) G-^^G-^ ($75 + $7(7+) G"^ ^ 
G^i ($7A) G^i ($79 + $79+) G^^ ($79+) + G^^ ($79) G^^ ($7A + $7A+) G^^ ($79+) 
+G-1 ($79) G^i ($79 + $79+) G-^ ($7A+) 

4 G^i ($79) G-^ ($79) G^i ($79) G^i ($79+) + 4 G^i ($79+) G^^ ($79+) G^^ ($79+) G^^ ($79) 
+3 G-^ ($79+) G^i ($79+) G^i ($79) G^i ($79) 
+3 G-^ ($79+) G^i ($79) G-^ ($79+) G^i ($79) . 

(4.32d) 



Here A is given by Eq. (4.29), which we restate for completeness: 

f r G-U) 



A 



-2i{4G-^ 
-2i 



(4.32e) 



($79) G-i ($79+) + ($79+) ($79) + ($79+) ($7"?+) 

I 

For the effective action to 0(9"*), A is needed only to This action needs to be augmented by a prescription 
0(9^). for the role of ^. From the above discussion it follows that 

contractions of ^ with 9 are given by the 9-propagator as 
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it follows from A 



(2) 
off ' 



(1^912 S934) 



(^912 °4> 



is a causal function of the complex frequency z. G is a 
coupling constant whose bare value is 
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27T 

<5l3<524 ^12 H ;r^Si-2,3-4 — 

V 1 



'P12 fSi 

^2^2 



G = 1/wF 



(4.36c) 



while contractions of ^ with ^ are given by 



{J12 J34) - ^ ^1-2.3-4 ^ ^ 



Here we have defined 



Al-2 



-E 



(^1-2, 3-4 ¥'34 



(4.33a) 



(4.33b) 



(4.34) 



It constitutes the natural coupling constant for a loop 
expansion. Sd-i = 27r''/Vr(d/2) with r(a;) the Gamma 
function is the surface area of the unit {d — l)-sphere, 
and dflp is the angular integration measure for the unit 
vector p for fixed unit vector fc in d dimension s. Fo r 
general dimensions d > 1 the latter takes the fornP^IlSl 



di] 



'^^^=fmJ-i (l-ry2)W2 ' (^-^-1>0), 

(4.37a) 

where e = d — 1 and 77 = cos 6 with 9 = <{p,k). The 
limiting expression for d — 1 is 



Physically, is the hydrodynamic part of the density 
susceptibility per spin in Hartree-Fock approximation, 
see the discussion below. 



dQp = 



Equations (4.32| through (4.341 completely specify the From Eqs. (4.36b), (4.371 we obtain 



df] [6{r] + 1) + 5(77 - 1)] , (d = 1) . 

(4.37b) 



effective theory for the purposes of a loop expansion for q- 
correlation functions.'^ They are the central formal result 
of the present paper. Their derivation makes is clear how 
to derive the effective action to any desired power in q. 



2. Gaussian propagators 

In order to make the theory suitable for actual cal- 
culations we need to explicitl y eva luate the correlation 
function x^^^ defined in Eq. (4.34). If we absorb the 



2Fi(l,l/2,d/2;l/z2) , 



(4.38) 



where 2F1 is Gauss's hypergeometric function. In d = 
1, 2, 3 this reduces to the familiar expressions for the hy- 
drodynamic part of the Lindhard function in these di- 
mensions: 

(PdMz) = -izjiX-z") , (4.39a) 
^d=2[.z) = sgn(Imz)/v/l - z2 



Hartree-Fock self energy into the chemical potential, as 
explained in Sec. 
reads explicitly 



= VVz + iVz-1 



IV B 1 the integral that defines x 



(0) 



^.=3(.) = ^ In 



(4.39b) 
(4.39c) 



1 



(4.35) 



The integral can be performed exactly. However, this is 
not necessary for most purposes, as the crucial hydrody- 
namic structure of x*-^^ is preserved in the well-known 
approximation, valid for small values of fc = |A;| and 
|fl„|, that performs the radial part of the momentum in- 
tegral by means of a contour integration over the interval 
—00 < < cxD, which makes the frequency sum trivial."^ 
Within this scheme, which we refer to as the AGD ap- 
proximation, we are left with an angular integral only. 



X^°Hk,inr,) - -Nf ipa{Gin„/k) , (4.36a) 



The Gaussian propagators, Eqs. (4.33 1, are now explic- 
itly specified. 



3. The Goldstone propagator and the density susceptibility 

From the preceding two subsections we see that the 
Goldstone modes are given by the propagator ipii{z), Eq. 
(4.38), which determines the density susceptibility in 
saddle-point (i.e., Hartree-Fock) approximation via Eq. 
(4.36a). This is just the hydrodynamic part of the famil- 
iar Lindhard function, see Eqs. (14.39 1. It is illustrative 



to also consider the density susceptibility in the Gaus- 
sian approximation to the effective field theory. From Eq. 
(2.18a) we see that the density susceptibility x is given 
by a {q q) correlation function plus terms that result from 
the coupling of q to P. Eliminating that coupling by the 
shift that was discussed in Sec. IIV C II leads to Fermi- 
liquid corrections. Using Eq. (4.33a) we obtain for the 



where 



density susceptibility in Gaussian approximation 



—if 1 

ipd{z) = - — / dnp — 

J p ■ k z 



(4.36b) 



X(fc,i0„) 



1 



2x<°Hfc,^»«) 

(i + 27VFr)2 i + 2jx^^){k,inr. 



(4.40) 
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We recognize this as the random-phase approximation 
(RPA) for the density susceptibihtj*^ with additional 
Fermi-hquid corrections. One of its characteristic fea- 
tures is the collective mode known as zero sound that 
results from a real root of the denominator. To see this 
we define another causal function 



1 



(4.41) 



X can be expressed in terms of this as 

-2iVF Gfln 



{l + 2NFTy k 



ifi'-^'^Gin^/k) . (4.42) 



From Eq. (4.36a I we see that, apart from the Fermi-liquid 
corrections, (/s^ relates to ipd the same way x relates 
to x^°^- Using Eqs. (|4.38l |4.39l) in Eq. kAlh one finds 



(s) . 

that in d = 1 the spectrum of (p^ is exhausted by the 
zero-sound poles. In d > 1 a particle-hole continuum 
emerges, but the zero-sound poles outside of the contin- 
uum remain. This is illustrated in Fig. [T] Note that the 
zero-sound modes are collective density fluctuations that 
are the result of particle-number conservation, whereas 
the Goldstone modes result from the spontaneously bro- 
ken continuous symmetry discussed in Sec. |III| and are 
not related to a conservation law. This manifests itself, 
for instance, in the fact that the Goldstone modes acquire 
a mass at nonzero temperatures, whereas the zero-sound 
mode does not. 



Equation ( 4.42 ) provides an important identity that al- 



lows to identify the thermodynamic density susceptibility 
dn/d^ in terms of the parameters of the field theory. For 
simplicity we consider the case d = 1. Using Eqs. (4.39a 
and (|4.41|) in (|4.42|) yields 



Xik^iftn) 



-2Nf 



1 + 2NfT (1 + 2iVFr)fc2/G2 + 

(d = 1) . (4.43a) 

The compressibility sum rule,'^ 



then yields 



dn/dfi = — x(A; — 0,«ri„) 

dn _ 2Nf 
dfl ~ 1 + 2NfT ■ 



(4.43b) 



(4.44) 



This is consistent with another argument, dn/dfj, is re- 
lated to th e zero- sound speed s by = {n/m) /{dn/dn)^ 
From Eq. (|4.43 a|) we see that {1 + 2NfV)/G'^ , which 



also yields Eq. (4.44 1 



While we have derived this relation from the Gaussian 
theory, it is expected to be an exact identity that remains 
true if dn/dfi and F are replaced by their renormalized 
counterparts. This is because the structure of the den- 
sity susceptibility, which is governed by particle-number 
conservation, must stay the same. This is an important 
point since it is known from perturbation theory that 
dn I d[i is not singularly renormalized near d = 1 




FIG. 1; The spectrum <(S^''"(a;) = V^&^''^\z = -I- iO) for 
2Np^ — 0.5 as a function of the frequency lj in dimensions 
d = 1, d = 1.5, and d = 3. The vertical lines denote 5-function 
contributions to the spectrum. 



E. Density formulation of the effective field theory 



As we have mentioned after Eq. (4.151 and (4.31), the 



current theory explicitly keeps all of the soft modes that 
were identified by the Ward identity in Sec. |III| and hence 
is local. This is in contrast to the theory formulated 
in Ref. 1T51 which was a formulation in terms of s-wave 
{i = 0) modes only. Technically this is achieved by con- 
straining Bn mjx, x) , whose spatial Fourier transform is 
given by Eq. (2.10a), to a local matrix Qnm{x) by means 
of a Lagrange multiplier field Anmix), and formulating 
the theory in terms of Q{x) and A(a;). This amounts 
to integrating out the modes in the higher angular mo- 
mentum channels, and since these are soft the result is a 
nonlocal field theory. While this is undesirable for many 



18 



purposes, for certain calculations it may be advantageous 
since it makes it easier to see which loop integrals as in- 
frared divergent and which are not. The nonlocality, on 
the other hand, makes no difference in practice unless one 
is explicitly interested in angular momentum channels 
with ^ > 1 or wishes to calculate non-local correlation 
functions. In the remainder of the section we therefore 
present the theory in such a density formulation. 

The Fermi-liquid saddle-point solution is the same as 
in Sec. |IV A 1[ except that the saddle-point value of Q 



is a position- independent object Qn given by Eq. (4.3a 



summed over the wave vector k. The further develop- 
ment proceeds as in the phase-space formulation, except 
that the various products of Green functions or their in- 
verses get replaced by products of convolutions of Green 
functions in wave-vector space and inverses of such con- 
volutions. The frequency structure of the theory remains 
unchanged. The Gaussian action, after decoupling P and 
q, then reads 



r=0,3 1.2 

3,4 



324 



24 



r^l2(^)r^34( 



E E^i4E [°^12W ('5l3'524^^-<5i_2,3-42Tr) 5!^34 (-fc) - ^13 ^524 ° A,2 A34 (-fc) 



r=0,3 1.2 

3,4 



(4.45) 



This replaces Eq. (4.25). In this subsection, 1 = rii, etc. represents Matusbara frequency labels only, in contrast to 



the notation used in the other parts of this paper. We have defined 

•^12 (fc) - :j^5]G„,(p)G„,(p-fc) 



$i2(fc) if (ni + l/2)(n2 + l/2) >0 
ifiuik) if (ni + l/2)(n2 + l/2) <0 



(4.46a) 



which replaces Eq. (4.12c I. Similarly, the products of inverse Green functions in the higher terms in the fluctuation 



expansion get replaced by more complex expressions, which involve convolutions of Green functions, with the same 
net power of Green functions. To account for this, we define 



$J'"|^j(fei, . . . , k„i) if ni through rim all have the same sign 
¥'i™m('''i! ■ • • I ^m) if through rijn do not all have the same sign , 



(4.46b) 



The cubic and quartic terms, instead of Eqs. (4.18 4.19), now take the form 

-a 



y2 



E E ^'=1 

1,2,3 k1.k2.k3 



+k2+k3,0 



x[^f^,{ki,k2) $r2'(fel)1>2-3'(fc2) $3"l'(fe3) {tr [(P - A) 12 (fcl ) (P + P+ - A - A+)23(fe2)(P+ " A+)3l(fc3)] 

+^Skfci> ^2) $r2'(fei) 'P2iik2) ^31(^3) tr [(P + P+ - A - A+)i2(fci) (^23(^2)^^1(^3) + ^^3(^2)^31(^3))] } • 

(4.47a) 



(4.47b) 
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with 



1,2,3,4 kiMMM 

X [4tr [{P - A)i2(fci)(P - A)23(fe2)(P - A)34(fc3)(P - A)|i(fc4)] 

+4tr [(P - A)+ (fci)(^ - A)+3(fc2)(P - A)+ (fc3)(P - A)4i(fe4)] 
+3tr [(P - A)i2(fci)(P - A)23(fc2)(P - A)+ (fc3)(P - A)+ (fc4)] 

+3tr [{P - A)i2(fcl)(P - A)+3(fc2)(P - A)34(fc3)(P - A)+ (fe4)] 



(4.47c) 



X! '5fei+fe2+fc3+fc4,0 l^ltLC^l' ^^2, ^3 

1,2,3,4 fcl,fc2>3>4 



: {$r2'(fcl) *2V(fc2) ^34'(fc3) ¥'ri'(fc4) tr [{P + P+ - A - A+)i2(fcl)(P + P+ - A - A+)23(fc2) 

X (^34(^3) tiliki) + ^tii^s) ^4l(fc4))] 
+ $r2'(fcl) ^2iik2) <f34'(fe3) ¥'ri'(fc4) tr [(P + P+ - A - A+)i2(fcl) ^(fcs) 

X (P + P+ - A - A+)34(fc3) ^li(fc4)] } , (4.47d) 



^-4^°'''^ ^ ^ H H 4i+fc2+fc,,+fc4,0¥'m4(^l'fc2,fc3)'Pl2^(fcl)'/'23^(fe2)'y534^(fc3) 9541^(^4) 

1,2,3,4 fci,fc2)'S3i'S4 

xtr [^12(^1)^^3(^2)^34(^3)^11(^4)] ■ 



The relation between P and P is 

:?Pl2(fe) - 0Pi2(fc) - $i2(fc) 27T^ 5i_2,3-4 °g34(fe) : 



(4.47e) 



(4.48) 



3,4 



which replaces Eq. (4.24a I. Note that the vertices in Eqs. (4.47) are invariant under cyclic permutations of the 



frequency and wave vector indices. Also note that these vertices are in general not finite in the hydrodynamic limit 

of small wave vectors and frequencies, since the n-point susceptibility tp^^-' ni^i' ■ • • > ^n) diverges in this limit .E^ This 
formulation therefore does not yield a local field theory. 



P and A can now be eliminated in analogy to Sec. IV C 2 The identities (4.28) turn into 



(P-A)i2(fc) = $12(fc)(rP)i2(fc) 
= cI>i2(fc)(7A)i2(fc) , 



with 



(4.49a) 
(4.49b) 



(fP)i2(fe) = 2Tr^5i_2,3-4P34(fe) . 

3,4 

(7A) is defined analogously, and so is (75). A in terms of q takes the form 

Ai2(A;) = ^ ^ fillip, k-p) ip^^{p) ip^^{k - p) [ii^ip) ^^2!^ - P) + itzip) h2{k - p)] 

3 P 

-2* E ^ E *^32(P, - P) [(7'Z)i3(p) (7g+)32(fc -P) + (79+)l3(p)(79)32(fc - P) 

+ (7'?)i3(p)(79)32(fc - P)] + 0[q^) , (4.50) 



3 ^ P 



which replaces Eq. (4.291 



Note that here, and in Eqs. (4.51) below, 79 always arises from a P that was shifted by means of Eq. (4.48), and 



therefore carries the same implicit frequency restrictions as P and A. That is, (7g)i2(fe) implies (ni + l/2)(n2 + l/2) > 
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0. A related remark is that the notation <(5i23' ¥'1234' ^^c. just specifies that the convolution has hydrodynamic content, 

1. e., that all of the frequencies do not have the same sign, but not which frequencies are positive or negative. This 
information, which is necessary for calculating physical quantities, can be obtained by an inspection of the matrix 
fields that multiply the convolution. 



The effective action analogous to Eqs. (4.32) then reads 

Aff - 4ff^ + AAI^ + a4^' + 0(g5) . (4.51a) 

with a Gaussian part 

4ff = -8 E E ^ E "9i2('^) - ^1-2,3-4 2r7) "gg^l-fc) (4.51b) 



r=0,3 1,2 k ^ -my ) 



and nonlinearities 



1,2 fci,fc2 



3 (=3 

r ( 

X 



{f'vl-ii^iM) ¥'23^(^2) V53i^(fc3)tr [(79 + 7'7"^)i2(fei) (fo(fc2) ^^1(^3) + ^^3(^2) ^31(^3))] 
+$f23(fci>fe2)tr [(79)12(^1) (7g + 79+)23(fc2)(79+)3i(fc3)]} , (4.51c) 



-^iff = ^^^E E '^fei+fc2+fe3+fc4,0 </'lt34(^l'fc2,fe3)V'l2^(fcl)¥'23Hfc2) 9534^^3) <^«(fc4) 



1/3 

1,2 fci,fc2 
3,4 fc3,fc4 



4i 



Xtr [^12(fei) ^J3(fe2) ^34(^3) ^ii(fc4)] 



1/2 

1,2 fci ,fe. 



E E '^fci+fc2+fc3,0<Pm('^l''^2)</323nfe2)<y53i^(fc3) 

Xtr [(7A + 7A+)l2(fcl) (fo(fc2) ^i(fe3) + ^+(fc2) ^3l(fe3))] 

2E^Etr [Ai2(fc)7Aji(fc)] 



12 ^ fe 

E E ^fel+fe2+fc3+fc4,0 '/^itLC^l' ^^2, ^3) 

1,2 fci,fc2 



3,4 k^Mi 
X 



{</'34^(fe3) 9541^(^4) tr [(79 + 7g+)i2(fci) (79 + 7'7+)23(fc2) (^34(^3)^41(^4) + ^^4(^3) ^4l(fc4))] 
+V'^3^('S2) <^r/(fc4) tr [(79 + 7g+)i2(fci) ^23(^2) (79 + 79^)34(^3) it\K^^)\ } 
E E ^fei+fc2+fc3,o$l2U'^i''^2,fc3){tr [(7A)i2(fci)(7a + 7'Z+)23(fc2)(79+)3i(fc3)] 



4i 

1,2 fci,fc9 
3 fc3 - 



+tl- [(79)l2(fcl) (7A + 7A+)23(fc2) (79+)3l(fc3)] 
+tr [(7g)i2(fci) (79 + 7g+)23(fc2) (7A+)3l(fc3)] } 

E E '^fei+fc2+fe3+/=4,o*m4('«i'^2,fe3) {4tr [(7'7)i2(fci) (7'7)23(fc2 (79)34(^3) (79^)41(^4)] 



1/3 

1,2 fci,fe2 
3,4 fc3,fc4 



+4tr [(7g+)i2(fei) (7'Z+)23(fc2 (7'7+)34(fe3) (79)41(^4)] 
+3tr [(79)12(^1) (79)23(^2 (7g+)34(fc3) (7g+)4i(fe4)] 

+3tr [(7g)i2(fci) (7g+)23(fe2 (79)34(^3) (7g+)4i(fe4)] } , (4.51d) 



where A in terms of q is given by Eq. (4.50 ). The rules about treating ^ are the same as in the phase space formulation. 
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and the Gaussian contractions are 



V 
16 



^13 '524¥'l2(fc) + 27T(5i_2,3-4 



^ujk) y34(fc) 

l + 27xi%(fc)J 



yi2(fc) V?34(fc) 

l + 27xi%(fe) 



where 



(;l°22(fc) = -T^(5i_2.3-4^34(fe) . 



(4.52a) 
(4.52b) 

(4.52c) 
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The above expressions provide aU of the information 
necessary for one-loop renormahzations of two-point cor- 
relation functions in the density formalism. 



V. A SIMPLE APPLICATION: PERTURBATION 
THEORY FOR THE DENSITY OF STATES 

As a simple demonstration of how the formalism de- 
veloped above works, let us consider the density of states 
(DOS ). We first use th e phase-space formalism from Sec. 
IV Dl From Eq. (12. 171) we see that the DOS is given by 



N{uj) 



Np + ReQ{iuj„ 
Nf + SN{uj) , 



iO) 



(5.1a) 



where Np is the DOS at the Fermi level including the 
interaction in Hartrec-Fock approximation, and 



4 1 



(5.1b) 



where we neglect 2-loop contributions and less divergent 
terms. Inserting the Gaussian propagator, Eq. (4.331, 
and choosing ni > 0, we find 



-2ij 1 

TT V 



m<0 



l + 27xi°-„.(fe) 



(5.4) 



Here (^'■^)++'~ is given by Eq. (4.46b I, and we have ex- 
plicitly indicated that the first two frequency arguments 
are positive while the third one is negative, and x*^*^' is 
given by Eq. (4.52c I. The same result is obtained if we 
realize that 6N is given by {Pnn{p)) in the density for- 
malism of Sec. IV E and use Eqs. (4.49a I, (4.501, and 
(|4!52bl). 



To find the leading hydrodynamic contribution to Q 
or SN we can calculate v?^'^^ in the AGD approximation, 
whence it depends only on ni — m. Performing the in- 
tegral atr = a; = Oind=l + e dimensions, we obtain 



We next observe, from Eq. (4.28) 



(A 



11/ 



22/ 



(5.2) 



The second contribution to (P) is proportional to the 
first one with an additional frequency integration which 
decreases the effects of the soft modes. The leading hy- 
drodynamic contribution to Q is therefore given by Eq. 
(5.1b) with P replaced by A. Using Eq. (4.29) we can 



express Q in terms of a loop expansion in terms of the 
g-correlation functions. To one-loop order the expression 
for A in Eq. (4.29) suffices. Of the two contributions to 
(All), the second one has no hydrodynamic content due 
to the frequency restrictions inherent in $ or 75. From 
the first contribution we obtain 



ttV 



(5.3) 



SNiLo = 0)/Np = — [1 + 0(6)] ^ f{Al) + 0{G^) , 
e TT 

(5.5a) 

where G is the loop expansion parameter from Eq. 
(|436c|. 



and 



Al = 2iVF7 = 



2NfT 
1 + 2NFr 



(5.5b) 



(5.5c) 



pi 



is the singlet s-wave Landau scattering amplitude in the 
particle-hole channel. 

There are several interesting aspects of this result. The 
1 /e singularity reflects the instabilit y of t he Fermi-liquid 
against a Luttinger liquid in d = j^lMllHII] ^j^g current 
formalism emphasizes that this is a result of the Gold- 
stone fluctuations getting so strong that they destroy the 
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ordered phase. In fixed dimension 1 < d < 3 as a func- 
tion of the imaginary frequency the corresponding result 
is 



(Q(zw„) - Q(zO))/iVF = cd /d(A^o)< 



d-l 



(5.6a) 



where > interpolates smoothly between cx — 
1) for d — > 1 and Cd oc 1/(3— d) for d — > 3. /d is a function 



that continuously evolves to / given in Eq. ( 5.5b I as d — >■ 



1, is positive definite for 1 < d < 3, and has the property 
fd[x — >■ 0) oc x^. For d = 3 one finds 



(g(zw„) - Qm)/NY 



23 /3(^o)^nlog(Vw„) 

+0{^l) (5.6b) 



with C3 > 0. 

After analytically continuing to the real axis these re- 
sults imply for the leading nonanalytic frequency correc- 
tion to the density of states 

5N{u)INf ^nd{2-d)\uj\'^'^ (5.7a) 

for 1 < d < 3, and 

5N{lo)/Nf = -n^uj'^\og{l/\uj\) (5.7b) 

for d = 3. Ud is positive definite and interpolates 
smoothly between Ud oc l/(d — 1) for d — >■ 1 and 
rid (x. 1/(3 — d) for d — >■ 3, and > 0. 

Note that to one- loop order there is no |a;| contribution 
to the DOS in d = 2, in agreement wit h earl ier results ob- 
tained with diagrammatic techniques .^^SlSMIll Equations 
(5.5a), (5.6 1, and (5.7) are analogous to the C oulomb 
or zero-bias anomaly in a disordered Fermi liquid .fSIlSIini 
There, the singularity at zero frequency is proportional 
to l/(d — 2), and the nonanalytic frequency dependence 
is The exponents of the respective frequency 

nonanalyticities reflect the scale dimensions of the least 
irrelevant operator at the disordered and clean Fermi- 
liquid fixed point, which are — (d — 2) and — (d — 1), and 
the corresponding dyna mical ex ponents, which are z = 2 
and z = 1, respectively.l^SlIlIIIl If -^ve denote the least ir- 
relevant operator by u, this leads to a homogeneity law 
in the clean case 



6N{uj) = 5iV(a;6,u6~('^"i') , 



(5.8) 



which immediately leads to the frequency dependence re- 
flected in Eq. (5.7a I. The sign of the effect at zero fre- 
quency, Eq. (5.5a), reflects the fact that a repulsive in- 



teraction will lead to a decrease of the density of states 
at the Fermi level. 

The dependence of the anomaly on the interaction am- 
plitude also warrants a comment. The function /(x), 
Eq. (5.5b I, is positive definite for < a; < 1, and 



f{x — 0) oc x^. The effects of the soft modes thus ap- 
pear only at second order in the interaction, whereas Eq. 
(5.3) or (5.4) naively suggest an effect at first order. The 



technical reason is that the integral in Eq. (5.4) vanishes 
if one puts 7 = in the denominator of the integrand. 



(a) 



(b) 



(c) 



FIG. 2: Diagrammatic contributions to the DOS at (a) first 
and (b) second order in the interaction amplitude, (c) Shows 
a ladder resummation of the diagram in (a) plus the first 
diagram in (b), which is identical with Eq. (|5.4[). 



This is consistent with many-body perturbation theory, 
as it must be: The diagram in Fig. [2] (a), which is the 
only contribution to the DOS to first order in the inter- 
action (Hartree diagrams have been absorbed into the 
Green function) , does not allow for a mixing of retarded 
and advanced degrees of freedom and hence has no hy- 
drodynamic content. The first diagram in Fig.[2|(b) does, 
and its hydrodynamic part is identical with Eq. ( |5.4 ) to 
second order in 7. The hydrodynamic part of its lad- 
der resummation. Fig. [2] (c), is identical with the full 
Eq. (5.4). The other contribution to second order in the 
interaction, i.e., the second diagram in Fig. [2] (b), rep- 
resents a particle-particle channel contribution. In the 
present formalism it will appear if a p articl e-pa rticle in- 
teraction amplitude is added to Eqs. (2.4b) or (2.13).^ 



Note that the effective theory picks out only the hydrody- 
namic parts of these diagrams, whereas the many-body 
diagrams contain non-hydrodynamic pieces as well. Also 
note that the current theory, at any given order in the 
loop expansion, in valid to all orders in the interaction 
amplitude, whereas many-body perturbation theory by 
necessity involves an expansion in powers of the interac- 
tion. 

In addition to the above discussion there are other in- 
teresting aspects of the DOS that include the behavior 
in the case of a long-ranged Coulomb interaction,'^^ Sll 
as well as the two-point correlation function of the local 
DOS. These effects will be discussed elsewhere.'^ 



VI. SUMMARY, DISCUSSION, AND OUTLOOK 

A. Summary and Discussion 

We have derived an effective field theory that de- 
scribes the soft modes, and their effects on observables, 
in electron systems without quenched disorder. The soft 
modes have been identified by means of a Ward idendity; 
they are particle-hole excitations with a linear frequency- 
momentum relation. They are the Goldstone modes of 



23 



a continuous symmetry, namely, a rotation in Matsub- 
ara frequency space that relates retarded and advanced 
degrees of freedom and is spontaneously broken when- 
ever the quasiparticle spectral weight is nonzero. This 
identifies a Fermi liquid as an ordered phase with the 
quasiparticle spectral weight as the order parameter and 
the soft particle-hole excitations the corresponding Gold- 
stone modes. A nonzero temperature the Goldstone 
modes acquire a mass, and so do subsets of them in the 
presence of a magnetic field or other symmetry-breaking 
external fields. Quenched disorder restricts the Gold- 
stone modes to the s-wave or £ = channel, whereas 
all higher angular momentum channels become massive. 
This is the reason why the clean theory is more com- 
plicated than the corresponding one for disordered elec- 
trons: There are many more soft modes to take into ac- 
count in clean systems. It is important to remember 
that the Goldstone modes are not density excitations, 
and their existence is not related to particle-number con- 
servation; the density susceptibility is characterized by 
zero-sound excitations regardless of whether the symme- 
try is broken or not, and at any temperature. 

Technically, we have borrowed heavily from techniques 
and concepts developed by Wegnei^^ for disordered non- 
interacting electrons and generalized by Finkelsteirfi2] to 
the interacting disordered case. However, the resulting 
effective field theory in the clean case does not take the 
form of a generalized nonlinear sigma model, as it does 
in the disordered case, but is significantly more complex. 
This complexity reflects the presence of many more soft 
modes. One manifestation of this complexity is the fact 
that the effective action contains terms to higher than lin- 
ear order in the interaction amplitude, whereas in the dis- 
ordered case only linear terms appear. The reason is that 
in the clean case the Green function is soft, and the ad- 
ditional powers of the frequency that accompany higher 
powers of the interaction are compensated by multiplica- 
tive Green functions. In the disordered case, in contrast, 
the Green function is massive because of the elastic re- 
laxation rate, and terms of higher than linear order in 
the interaction are irrelevant in a renormalization-group 
sense. 

We hav e follow ed Wegneif^ and later work based on 
his ideaJsHnmilll bosonizing the theory by introduc- 
ing classical matrix fields that are isomorphic to bilin- 
ear products of fermion fields. Technically, this is done 
by introducing a Lagrange multiplier constraint and in- 
tegrating out the fermions. This approach is very differ- 
ent from that of Shankar,'23who applied renormalization- 
group techniques directly to the fermionic theory to de- 
rive Landau Fermi-liquid theory. The Lagrange multi- 
plier field, whose expectation value plays the role of a 
self energy, is soft and needs to be kept to all orders to 
avoid introducing spurious soft modes. This has been 
achieved in terms of diagram rules for the effective the- 
ory, with the role of the Lagrange multiplier being to 
eliminate the noninteracting part of certain propagators. 
The effective action has not been given in closed form. 



but an explicit method has been devised to construct 
it to any desired power in the fundamental soft-mode 
field. This allows for a systematic loop expansion, with 
the inverse Fermi velocity playing the role of the bare 
loop expansion parameter. The interaction, by contrast, 
is not treated perturbatively. Order by order in a loop 
expansion, the theory thus produces the leading hydro- 
dynamic (i.e., long- wavelength and low-frequency) effects 
to all orders in the interaction, and the interaction is not 
required to be weak in any sense. This, together with 
its ability to isolate the hydrodynamic contributions to 
any desired observable, sets the current theory apart from 
many-body perturbation theory which expands in powers 
of the interaction, with any resummations to all orders in 
the interaction being opportunistic and not controlled. 

A related feature of the current theory that is even 
more important is its amenability to renormalization- 
group techniques. Since it is formulated in terms of 
an infinite power series in the fundamental field, with 
each term characterized by an n-point vertex and a cor- 
responding correlation function, the structure of which 
is dictated by symmetry, the bare theory also suggests 
the structure of the corresponding renormalized theory, 
which many-body perturbation theory cannot do. This 
feature has been exploited extensively in the disordered 
case to describe, i nter alia , the Anderson-Mott metal- 
insulator transition,'i^lIlllII which had proven very hard 
to do by means of many-body theory. 

Certain aspects of the present theory had already been 
present in the theory put forward in Ref. il6t However, 
since the latter focused on disordered systems, it was for- 
mulated in terms of s-wave variables, with the higher an- 
gular momentum modes effectively being integrated out. 
In an application to the clean case, this led to a nonlocal 
theory, i.e., the 3-point and higher vertices did not exist 
in the hydrodynamic limit of small frequencies and wave 
numbers. In addition, no complete separation of soft 
and massive modes was achieved in that theory for the 
clean case. Both of these points have been remedied in 
the present formulation. The issue of the higher-angular- 
momentum modes is handled by formulating the theory 
in terms of phase-space variables, which explicitly keeps 
all angular momentum channels, and the separation of 
soft and massive modes has been achieved by expressing 
the massive modes in terms of the soft ones in a power- 
series expansion, albeit not in a closed form, in contrast 
to the disordered case. The relation to the formulation in 
terms of s-wave variables, which is still useful for certain 
calculational purposes, has been explained in Sec. IIVE] 



B. Outlook 

Formalism for formalism's sake is not very useful, and 
a crucial test for the present theory will be any physical 
problems it can solve more easily than other methods. 
We conclude this paper by listing several applications 
that we plan to pursue in the near future, many of which 
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we already have obtained preliminary results for. 

The density-of-states calculation in Sec. |V] is meant 
only as an illustration of how the loop-expansion-based 
perturbation theory works within the present context; it 
does not achieve anything that cannot easily be done 
with many-body diagrams. It nevertheless has some 
important advantages, and points to some crucial fu- 
ture developments. For instance, as was shown in Ref. 
rroi the Fermi-liquid state is represented by a stable 
renormalization-group fixed point, and the leading irrel- 
evant operators with respect to that fixed point can be 
easily identified. This shows that the nonanalytic en- 
ergy or frequency dependence of the DOS near the Fermi 
level, Eqs. (5.6 1, is qualitatively exact; there cannot be 



ena, i.e., order in the spin-triplet channel. In Sec. IV A 2 



we have constructed a saddle-point solution that repro- 
duces the Stoner theory of ferromagnetism. By expand- 
ing about that saddle point, analogously to the expansion 
about the Fermi-liquid saddle point in Sec. IV B one can 



any contributions with a lower power of the frequency. 
This would be impossible to establish within many-body 
theory. A related point is that the field theory, in contrast 
to many-body perturbation theory, is not perturbative in 
the interaction; the results at any order in the loop ex- 
pansion are valid to all orders in the interaction. Other 
observables will also display universal nonanalytic behav- 
ior in the Fermi-liquid phase that can be understood as 
corrections to scaling at the Fermi- liquid fixed point, and 
they can be identified by the same techniques. One in- 
teresting quantity is the two-point correlation function 
of the local DOS. Since the DOS is the (zeroth moment 
of) the order parameter for the Fermi liquid, this plays 
the role of the order parameter susceptibility, and it is 
expected to be long-ranged, i.e., to exhibit a universal 
divergence for small wave numbers, everywhere in the 
Fermi-liquid phase. This expected phenomenon is analo- 
gous to a well-known feature of the magnetic susceptibil- 
ity in the ordered phase of a Heisenberg ferromagnet.'-^ In 
Ref. ISniwe have used scaling arguments combined with a 
preliminary version of the current theory to predict that 
the DOS-DOS correlation function diverges as T/fc^^"^ 
for small wave numbers k. We will revisit this problem, 
together with other aspects of the DOS, in more detail 
shortly. 

In addition, the effective theory suggests the exis- 
tance of a critical fixed point that describes a transition 
from a Fermi liquid to a non-Fermi-liquid state. This 
will be analogous in many respects to the critical fixed 
points that describe various versions of the Anderson- 
Mott transition.lSl The basic physical idea is that the 
fiuctuations in a Fermi liquid that are represented by 
the Goldstone modes become stronger with decreasing 
dimensionality and increasing interaction strength, and 
ultimately give rise to an instability of the Fermi liquid 
since the system can lower its free energy by going into a 
state that restores the symmetry and thus eliminates the 
Goldstone modes. A scaling theory for such a transition 
has been developed in Ref. (HUl and arguments have been 
given that suggest that the corresponding fixed point will 
appear in the present theory at 2-loop order. The loop 
expansion will then allow for a controlled description of 
the transition in d = 1 + e dimensions. 

The theory can also be applied to magnetic phenom- 



construct a systematic theory for quantum Heisenberg 
ferromagnets that improves upon and replaces the Hertz- 
Millis theory.l^^Illl Such a theory has been developed in 
Ref. ESI and it has been shown to have very interesting 
consequences, including the fact that the phase transi- 
tion at sufficiently low temperatures in a Heisenberg fer- 
romagnet is generically of first order, in agreement with 
experimental observations. This theory was based on a 
phenomenological treatment of the fermionic degrees of 
freedom combined with symmetry arguments, and it does 
not have the correct limiting behavior as c? — ^ 1. The cur- 
rent theory will allow for a systematic derivation valid in 
any dimension d > 1 that will serve as an important 
check and will probably point to new developments. It 
can also be applied to ferrimagnets, the quantum aspects 
of which have not been considered so far. Finally, if one 
considers band electrons instead of the nearly-free elec- 
tron model we have employed for simplicity, the theory 
can also describe antiferromagnetism. 

Treatments of the spin-triplet channel are not re- 
stricted to s-wave order. By keeping a triplet interaction 
amplitude in the p-wave channel, see Eq. (2.13d), one can 
describe a magnetic nematic or p-wave ferromagnet. In 
Ref. [2ni we have shown that the quantum phase transi- 
tion from a Fermi liquid to such magnetic state is gener- 
ically of first order, for reasons that are similar to those 
for the analogous statement for the s-wave ferromagnetic 
transition. This is in contrast to Hertz-type theorie^^ 
that predict a continuous transition. This result was also 
based on a preliminary version of the current theory, and 
the present complete formulation will allow for more de- 
tailed studies of exotic magnetic states. The same is true 
for superconducting states, in either the s-wave or higher 
angular momentum channels. 

All of the above examples involve phenomena that cru- 
cially rely on the coupling of some order parameter to 
fermionic soft modes. The present effective theory pro- 
vides a general framework for describing the universal 
aspects of ordered phases, and the quantum phase tran- 
sitions between them, of clean electron systems. 
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Appendix A: Effects of quenched disorder 

In this appendix we show that quenched disorder mod- 
ifies the Ward identity we discussed in Sec. |III| such that 
only density or £ = modes remain soft. 

Consider a static random potential u{x) that couples 
to the electron number density and has a Gaussian dis- 
tribution with a second moment 



{u{x) ^U{x~ y)/TTNF 



(Al) 



where {. . .}^^^ denotes the disorder average and U is di- 
mensionally an energy density. We then need to add to 
Eq. ( |2.13a ) a contribution 



-4dis = ^5]V(q)^tr(Q(fc;q))tr(Q(p;-q)) 

^ ^ q k 



~T'T 

q k,p 



C/(fc-p)tr {Q{k-q)Q{p--q)) . 

(A2) 

The trace now includes tracing over replica indices to 
handle the disorder average. Now we anticipate that the 
center-of-mass momenta fc and p will be pinned to the 
Fermi surface, and expand U{k~p) in Legendre polyno- 
mials. We then can write the disorder part of the action 



1^11 

P=0 T^rol ^ m= 



xtr (Q^'"(g)Q^-'"(-q)) . (A3a) 

Here we have expanded the Q-niatrices in spherical har- 
monics Y^™, 

Qi2\q) = E ^"("'^) q) , (A3b) 

fc 

and have defined relaxation rates 1/ri = U{q — 0) and 



1 

'rel 



27r 



d7jPi{r])U{kFV^{l~Tl)) (A3c) 



in the angular momentum channels. 

The derivatio n of t he Ward identity as in Sec. |III| now 
still yields Eq. ( |3.12[ ) (note Q° ° = Q'-°^), but it gives no 
information about the angular momentum channels with 
£ > 1. This is consistent with the results of perturbation 
theory. The Fermi-liquid saddle-point solution from Sec. 
|IV A 1| needs to be augmented by a disorder-induced self 
energy and becomes identical with the saddle point con- 
sidered in Ref. [161 An expansion about the saddle point 
to Gaussian order yields 



W!f(q = 0)Qlf(q 



l/^(fco)_^, (£) 



(A4) 



In order for the theory to be stable we need to require 
l/rj^^[~°^ > 1/t^^o1^°\ which will generically be the case. 
If it is not one needs to choose a saddle point in the 
angular-momentum channel that has the largest value of 
1 /t^^ci"^'' . We see that all correlation functions other than 
the one for £ = are massive, which is consistent with 
the information obtained from the Ward identity. 



Appendix B: An illustrative example: Soft modes in 
the classical 0(2) vector model 

In this appendix we discuss a classical 0(2) vector 
model as an illustrative example of how to separate soft 
and massive modes. Consider a 2-component classical 
field <p{x) = {Tr{x), a{x)) and an action 



dx 



lc^\x)+'-iVc^ix)f + '^{<j>'ix)f 

(Bl) 

with (V0)^ = di(f>jd^<j>^ . The partition function is 

Z = J D[(P] e--^['^l . (B2) 

The saddle-point equations for homogeneous field config- 
urations (tt, a) read 



rvr -|- UTT^ + UTTO^ 







(B3a) 

ra + ua^ + ucttt^ = . (B3b) 

For r < 0, a solution that minimizes the free energy is 

TTsp = , (B4a) 

CTsp = \J -rju = 00 • (B4b) 

Alternatively, we can choose TTgp to have any value with 
ItTspI < (/)o. Then 

a% = <g{\-^%l<g) (B5) 

also solves the saddle-point equations. This reflects the 
0(2) rotational symmetry of the model. 

Now consider fluctuations about the saddle point, a — 
(/io(l + (^cr), and redefine tt — (/>o7r, so that 8a and tt are 
both dimensionless. Then the action takes the form 



(B6a) 



where 



hi J dx {Vn f 



00 / dx tt'^ 



(B6b) 



AsaM = u^l J dx iSay + ^(f>l J dx {VSaY 

(B6c) 

-4i„t[7r, (5cr] = U0Q J dx 6a + j dx 7r'^{6a)'^ . 



(B6d) 
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(a) 



(b) 



int 



+ 



FIG. 3: Graphic representation of the contributions Asa and 
.4int to the action. Dashed and sohd lines represent the fields 
5a and vr, respectively. 



Here da and tt are functions of a;. A graphic representa- 
tion of the terms that depend on da is shown in Fig. [3j 

It is well known that the transverse fluctuation, i.e. the 
field TT, is the soft mode in this model.l^ tt is the Gold- 
stone mode that reflects the spontaneously broken 0{2) 
invariance, and it is governed by a Ward identity that 
ensures that the homogeneous transverse susceptibility 
diverges. The effective soft-mode theory is a nonlinear 
sigma model with an action 



A 



NL(tM 



dx 



and a partition function 



DM 



(B7b) 



Notice that the functional integration measure has 
changed by switching to tt as the only field. This is im- 
portant to preserve the symmetry of the problem. Tech- 
nically, terms arising from the measure cancel spurious 
mass terms that get generated in perturbation theory. 

In this simple case the nonlinear sigma model can be 
derived in closed form by writing ^ as a representation 
of 0{2) rotations, 



■k'^{x) -nix) 
—7r{x) ^1 — Tr'^{x) 





p{x) 



(B8) 



and integrating out the massive /9-fluctuations in saddle- 
point approximation (i.e., replacing p{x) by its saddle- 
point value (pQ, see Ref. [7j). For the more complicated 
matrix field theory we are interested in no closed-form 
derivation of the effective soft-mode theory has been 
found so far, but one can still derive the latter pertur- 
batively order by order in powers of the soft modes. To 
guide that derivation it is useful to consider the analo- 
gous process for the 0(2) model. In what follows, we 
thus pretend that the nonlinear sigma model were not 
known, and construct an effective soft-mode theory per- 
turbatively by considering Eqs. (B6) and systematically 
integrating out Sa. 



Returning to Eqs. (B6) we see that at the Gaussian 
level TT is indeed massless and Sa is massive. However, 



> 
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FIG. 4: (a) Tree-level and (b) one-loop contributions to the 
term of 0{'k^) in A^-^fyr] in an expansion in Sa loops, (c) 
Tree-level contribution to the term of 0{7v''). 



the 77^ term in Eq. (B6b| does not carry any gradients. 



If we were to simply neglect 6a this term would generate 
a mass for tt at the one-loop level. This mass violates 
the Ward identity and is spurious; it gets canceled if one 
properly integrates out Sa. The coupling between Sa 
and TT is thus crucial and needs to be kept. Formally 
integrating out Sa leads to an effective action in terms of 
TT only. 



where 



AA^[n] = - In J D[Sa] e--4.5-['5'^l--4-t[7r,5<T] 



(B9a) 



(B9b) 



We can now perform the Ja-integral perturbatively. 

Let us first integrate out Sa in a straightforward man- 
ner. The resulting loop expansion is an expansion in 
powers of l/cfiQ. The tree-level and one- loop diagrams 
that contribute to 0{-k^) are shown in Fig. [i] The tree 
diagram yields 



dx TT 



J dx (V7r2)2 (BIO) 



Combining Eq. (BIO) with Eq. (B6b) we see that the 



offending tt'' term cancels, and the gradient terms agree 
with the expansion of the nonlinear sigma model, Eq. 
(B7a), to order tt^. Terms of higher order in tt can be 



constructed analogously, see Fig. ^c) 

Alternatively, we can evaluate the integral in Eq. ( B9b ) 



by means of a saddle-point approximation for Sa for fixed 
TT. The saddle-point equation is shown graphically in Fig. 
[5] Analytically, we obtain 
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FIG. 5: Saddle-point equation for Sa for fixed tt. 



= ^7r^x)+(^l+^7T\x)^da,p{x) 



2 V'Saspix) + - iSa,^ix)r + - (5a,p(a;))^ . 

(Blla) 



In a gradient expansion, and keeping terms up to O(V^), 
this is solved by 



S(Tsp{x) = yj\ - 7r2(a;) - 1 



+ 



l-7r2(a;) 



vVi-^'(3;) + o(ff2v^ 



(Bllb) 



with 5 = c/2u(ly^. Note that higher powers of 5 are ac- 
companied by higher powers of grad ient s. Substituting 
this result back into the action, Eqs. ( B6 ) , we obtain the 



action of the nonlinear sigma model, Eq. (B7a 



In contrast to the (5cr-loop expansion above, this second 
method does not require knowledge of the 6a propaga- 
tor. It also involves only solving an algebraic equation, no 



integrals need to be performed. In contrast to the stan- 
dard derivation of the nonlinear sigma model, which uses 
Eq. ( B8 1 , it reqires no explicit knowledge of the sym- 
metry group and its representations. If no closed-form 
solution of the saddle-point equation could be found, one 
could still solve Eq. (Blla I iteratively, since it is inho- 
mogeneous. In this way one obviously can construct the 
effective action perturbatively, order by order in powers 
of vr. For an analysis of the effective action to any given 
order in a loop expansion this would be sufficient. This 
is the procedure we follow for the matrix field theory in 
Sec. HYCl 

Finally, we need to worry about the change of the in- 
tegration measure mentioned after Eqs. (B7). In the 



present formalism, this is generated by the Gaussian fluc- 
tuations about (5crsp. If we write 5a{x) = 5asp{x) -\- t{x) 
and integrate out t{x) in Gaussian approximation we 
obtain an additional contribution to the effective action. 
The complete effective action then is 



Acs = ^ I '^^ 



Tr \n^/l - tt'^ + 0{G), 



(B12) 



with G — ^/c(jy^ the coupling constant. This is the non- 
linear sigma model including the measure terms, see Eqs. 
dBTl). 
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